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INTRODUCTION 


The quadratic equation x 2 — x — 1 = 0 is known to have roots 

l±y/5 
X ~ 2 ’ 

but suppose that we decide instead to solve this equation using an iterative 
method. The equation can be rearranged as 

x = 1 + 

x 

so we could find successive approximations, x n , to a root by using the 
formula 

£n+1 = 1 "i • 

X n 

If we take x\ = 1 as our initial approximation, the sequence begins as 
follows. 


The roots of ax 2 + bx + c = 0 are 


x = 


—b ± V& 2 - 4ac 
2 a 


£1 
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+ 
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£1 
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+ 
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£2 

£4 

= 1 

+ 

1 




£3 


1 + 


1 _ 2 
T ~ T 

l 


3 

2 


1 + 


1 


1 + 


1 + 


3^5 — Id-— 1 + 

£4 


1 + 


1 + 


1 


*4 


It is interesting to see the Fibonacci numbers appearing as numerators and 


denominators of each x n , and we might conjecture that x n = 


F n +1 


, for all 


positive integers n. This is certainly true when n — 1, and if we assume that 
x k = then 


, , 1 , , Fk 

£fc+1 = H -— 1 + -= - 

Xk -Cfc+i 


Fk+ i + Fk Fk- (-2 


^fc+i 


Fk+ i 


So Mathematical Induction confirms that it is true that x n 
n > 1. 


F n 


+i 


, for all 


We met the Fibonacci numbers in 
Section 5 of Unit 2. 


The question remains as to whether or not the sequence x\, X 2 , £ 3 , ■ • • 
produced by this iterative process converges. The question is the same as 

F‘2 F 3 F4 

asking whether the sequence —, —, —, ... of ratios of successive 

F\ F 2 F 3 

Fibonacci numbers converges, and this problem has a nice geometric 
interpretation. Imagine squares with edge lengths 1, 1, 2, 3, 5, 8, 13, ... . 
These can be joined, in order, to form a sequence of rectangles where, at 
each stage, the new rectangle is formed by adjoining the next square to the 
longer edge of the existing rectangle. This process is illustrated in the 
following figure. Now the shape of each rectangle, as determined by the ratio 
of the longer to the shorter side, is the ratio of successive Fibonacci 
numbers. So the question becomes: does this sequence of rectangles have a 
limiting shape? 
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The answer, as we shall see from the work of this unit, is yes. If, for the 
moment, we assume that the sequence has a limit k, say, then it is not too 
difficult to see what k must be. For sufficiently large rectangles the ratio of 
the lengths of edges is very close to k, so we could take these lengths to be 
(approximately) ky and y. To obtain the next rectangle we add on a square 
of side ky as shown below. 


y + ky- 


y 

ky 

ky 

ky 


The ratio of the sides of this larger rectangle is also (approximately) k and so 

V + ky _ 
ky 


which, after cancelling the y, simplifies to k 2 - k - 1 = 0. Hence, if the 
sequence tends to a limit, it will be a root of the quadratic equation from 

which we started, x 2 - x - 1 = 0. The positive root is k = —, a 
number which is known as the golden ratio. 

The golden ratio was, along with 7r and y/2, one of the first irrational 
numbers known to, and used by, the ancient Greeks, arising originally in the 
geometrical problem of dividing a line segment AB at a point C such that 
AB AC , 

this common ratio turning out to be the golden ratio. 


The value of the golden ratio is 
approximately 1 . 6180 . 


A 

h 


C 

-+- 


B 

H 


The golden ratio crops up all over mathematics and is particularly 
prominent in problems concerning regular plane and solid figures. 
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As one example, it is the ratio of the radius of a circle to the edge-length of 
an inscribed regular decagon. The ancient Greeks appear to have been 
endeared with the golden ratio and it may have been deployed by their 
sculptors and architects for aesthetic purposes. It features in the structure of 
the Parthenon and the theatre at Epidavros; for example, the theatre has 21 
tiers of seats above the gangway and 34 below, 21 and 34 being consecutive 
Fibonacci numbers whose ratio is a good approximation to the golden ratio. 


Trigonometry shows that this ratio 


2sin(7r/10) ’ 

which is another expression for the 
golden ratio. 


Our main purpose in giving this example is to illustrate what is meant by a 
continued fraction. Each of the expressions 


1 + Id-j > 1 + 

1+ 1 


1 + 


1 + 


1 + 


1 + 


1 


1 + 


1 


1 

1 + T 


is called a finite continued fraction. Their limit is the infinite continued 
fraction: 

1 


1 + 


1 


1 + 


This is the simplest, and best-known example of a continued fraction. In the 
more general infinite continued fraction 
1 

ai H - j- 

0-2 d-j- 

a 3 d-j- 

a 4 H- 

the ai are positive integers, with the exception of the first one, a\, which 
may be a positive or negative integer or zero. 

The theory of continued fractions is extensive and their area of applicability 
very wide, so in this unit we can only hope to give you something of the 
flavour of them. We shall give a number of diverse applications and, in 
Unit 8, we shall see how they can be exploited to lead to a neat solution of a 
famous Diophantine equation. 

We shall first explore finite continued fractions and, in so doing, lay the 
foundations for our subsequent investigations of infinite continued fractions. 
We shall discover that every infinite continued fraction converges to an 
irrational number and, conversely, that every irrational number can be 
expressed, in a unique way, as an infinite continued fraction. 

We are all in the habit of writing down, and even manipulating, expressions 

such as 7 r, log e 3, etc., which represent irrational numbers. But 

when we come to practical calculations involving them we must deal with 
rational numbers. (Even today’s computers are limited in how many decimal 
places of 7r they can handle.) So the question of finding good rational 
approximations to irrational numbers is an important one. In the final 
section of the unit, we see how the continued fraction of an irrational 
number provides us with ‘best possible’ rational approximations. 
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1 FINITE CONTINUED FRACTIONS 


1.1 The continued fraction of a rational number 

We start with the key definitions. 


Definition 1.1 Finite Continued Fractions 


Let ai be any real number and let a 2 , a 3 ,...,a n be any positive real 
numbers. The expression 


1 

«i H- 

fl 2 H- 

a .3 H- 


1 

O n —l T 


is a finite continued fraction. The numbers ai, a 2 , ..., o„ are the 
partial quotients of the continued fraction. When all the partial 
quotients are integers the continued fraction is said to be simple. 


Continued fractions can be cumbersome objects and it is therefore 
advantageous to have a convenient symbol to display them. Any continued 
fraction is determined by its partial quotients and so we write the finite 
continued fraction in the definition as 


[ui, 02) O^) . . . ) On\. 


For example, 


[1,2,5] = ! + 



and the expression on the right simplifies to the rational number — which is 

the value of the continued fraction. Similarly the value of [-3,1,2,4] can be 
calculated as follows. 


[ — 3,1,2,4] = —3 + 


1 


1 + 



-3 + 


1 


>♦5 


= - 3+ u = 


30 

13 


Problem 1.1 _ 

Determine the rational numbers represented by the following finite 
continued fractions. 

(a) [2,3,4] (b) [0,2,2,2] 


Looking at the first step in the calculations in the examples above, note that 

|1,2 ' 5| = 1 + M “ d 3+pWj. 

In general, from the way the finite continued fraction is defined, we have the 
following property. 


Note that formally we think of the 
finite continued fraction and its 
value as being distinct. The 
continued fraction is an arithmetic 
representation of a number (its 
value). However we shall 
occasionally drop this distinction 
and write, for example, 
x = [—3,1,2,4] to mean that x is 
the value of this continued fraction. 
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First continued fraction identity 


[ai,02,03,... ,a„] — ai + 


[02, o3,..., a„ 


In our work we are going to be concerned almost exclusively with simple 
continued fractions; that is, ones in which all the partial quotients are 
integers. There will, however, be isolated occasions when it proves 
convenient to use continued fractions with real partial quotients for the 
purpose of providing general proofs, though the results will subsequently 
only be applied to cases involving simple continued fractions. 

In the text, for ease of notation, we shall use the abbreviation FCF to stand 
for a simple finite continued fraction. 

It is readily seen from the definition that the value of any FCF is a rational 
number. As it happens the converse also holds: any rational number can be 
expressed as an FCF. The Euclidean Algorithm provides a means of 
determining an FCF for a given rational number. We illustrate this method 
in the following example. 


Example 1.1 

225 

Determine an FCF for yyy. 

Dividing 225 by 157, the Euclidean Algorithm gives the sequence of 
equations on the left below. Each of these is written in an alternative form 
on the right. 



225 

, 68 

225 = 1 x 157 + 68 or 

— 

— 1 + 


157 

157 


157 

„ 21 

157 = 2 x 68 + 21 


— 2 + —— 


"68" 

68 


68 

„ 5 

68 = 3 x 21 + 5 


= 3 + — 


21 

21 


21 

1 

21 = 4x5+1 

Y 

= 4+ 5 


5 


5 = 5x1 


= 5 


1 



From the system of equations on the right, by continually using the next 
equation to substitute for the existing fraction, we obtain the following. 

225 


, 68 , 1 
157 + 157 + 157 

68 


1 + 


1 


„ 21 
2+ 68 


= 1 + 


1 


2+ J_ 

68 


21 


- 1 + 


= 1 + 


2 + 


2 + 


3 + 


21 


3 + 


21 

T 


= 1 + 


2 + 


3 + 


4 + 
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225 

In other words, = [1, 2,3,4,5], the partial quotients of the FCF being, in 
order, precisely the quotients as they arise in the Euclidean Algorithm. + 


Problem 1.2 _ 

Determine an FCF for each of the following rational numbers. 




This way of using the Euclidean Algorithm can be deployed to obtain an 
FCF for any rational number. Let us formalize this before trying more 
examples. 


Theorem 1.1 Equivalence of FCFs and rational numbers 

Every simple finite continued fraction has a rational value and, 
conversely, every rational number can be expressed as a simple finite 
continued fraction. 


Proof of Theorem 1.1 

We prove that the value of any FCF is rational by Mathematical Induction 
on the number of partial quotients. 

The induction is started by observing that any FCF with a single partial 
quotient, [di], is an integer and therefore rational. 

For the induction step, suppose that the value of any FCF with k partial 
quotients is rational and consider the FCF [d X , u 2 , a 3 ,..., a k+ 1 ] with k + 1 
partial quotients. We now use the first continued fraction identity: 

[di,d 2 ,d 3 ,...,dfc+i] = oi + 7 --- 7 . 

l«2> <*3, • • • > afc+ij 

Now ai, being an integer, is rational. Furthermore, by the induction 
hypothesis, [ 02 , <* 3 ) • • •, «fc+i] is rational and hence so is its reciprocal. It 
follows that [ai, 02 , d 3 ,..., dfc + i], being the sum of two rationals, is itself 
rational. 

Mathematical Induction therefore confirms that the value of any FCF is a 
rational number. 

Conversely, let - be any rational number. By the Euclidean Algorithm: 

a = qib + n, 0 <n<b or ^ = <Zi + ~r 

b b 

b = q 2 r\ + r 2 , 0 < r 2 < n or — = a 2 + — 

r 1 n 

n = g 3 r 2 + r 3 , 0 < r 3 < r 2 or — = q 3 + — 

r 2 r 2 


r n - 3 = Qn-i r n-2 + r„_ 1 , 0 < r n _i < r n _ 2 or = qn l + !iLJ. 

r n -2 r„_ 2 

r n -2 = q n r n -i +0 or -= q n 

Tn—l 


By definition [di] = di. 


As d 2 > 0, by definition, 

[d2, d 3 , . . . , dfc-j-l] ^ 0. 


When we write — for a rational 
0 

number we assume that d, b are 
integers with b positive and 
gcd(d, b) = 1 . 


The fact that the non-zero 
remainders decrease guarantees 
that the algorithm stops after a 
finite number of steps. 
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From the right-hand system of equations we obtain 


a 1 

b= qi + T =qi + 
r\ 


— 9l + 


92 + rT 
V2 


1 


<72 


1 

93 + Tjf 

^3 


= 9i + 


92 + 


93 + 


9n—1 H- 

Qn 


which shows that 7 has the FCF [ 9 i, 92 ,--, 9 n]- 
b 


Problem 1.3 - 

Obtain FCFs for the following rational numbers. 


/ s 65 , . 23 , . 

<*> 23 (b) 55 (C) 


-94 

13” 


The FCF of a rational number is not unique. In fact each rational number 
has exactly two FCFs. The ambiguity arises because of the freedom in 
choosing whether or not to terminate the FCF with partial quotient 1. To 
substantiate that remark, look again at Example 1.1 where we saw that 

225 1 

-= [1,2,3,4,5]. The final partial quotient, 5, can be written as 4 + -: 

157 t 

1 1 


+ 


4 + 


= ... + 


4 + 


1 


4 + 


1 


Consequently we have an alternative form for the FCF 
|P = [1,2,3,4,4,1]. 

Similarly, from Problem 1.3 part (a) above 
65 


23 


[2,1,4,1,3] = [2,1,4,1,2,1], 


the alternative being reached by writing the final partial quotient 3 as 2 + -. 
In general, if a n > 1, 

[<2l, (Z 2 ,. . • , On— 1 j <^n] = [<^1 ? <^2»• • • j 1 ) 1 ] ? 


since CLj « 


(^71 


1)+ r 


Problem 1.4 


Write each of the following as an FCF with an even number of partial 
quotients. 


(a) [1,2,3] (b) [-4,3,1,6,1] (c) [0] 


If an FCF has a final partial quotient of 1 then, in the way illustrated in 
Problem 1.4 part (b), the final 1 can be absorbed into the previous partial 
quotient. If we choose always to absorb such Is then the ambiguity in the 
FCF is removed and every rational number has a unique FCF whose final 
quotient is other than 1. (There is one exceptional number. The integer 1 
has FCF [1], where the final partial quotient 1 cannot be absorbed into the 
previous one because there is no previous one. The two FCFs of 1 are [1] and 
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[0,1], both of which end in 1.) We shall not give a proof of this uniqueness 
here for, in Section 3, we shall prove a general result concerning the 
uniqueness of infinite continued fractions which can be applied to the finite 
case. 


1.2 The convergents of a continued fraction 

Let us turn to the question of evaluating FCFs. First we need to introduce 
some further terminology. 


Notice that the final convergent C n 
is the value of the FCF itself. 

The reason for the name convergent 
As an illustration of the definition consider [1,2,4,3,2], It has five Wl11 become a PP aren t in due course, 

convergents. 


Definition 1.2 Convergents 

For each fc = l,2,...,nwe define the kth convergent, Ck, of the FCF 
[a\, a 2 ,... ,a n ] to be the value of [ai, 02 ,..., a*]. 


Cl = [ 1 ] = 1 

C 2 = [ 1.2]-1 + I»| 

C 3 — [ 1 , 2, 4] = 1 4 —y = 

2+- y 
4 

C 4 = [1,2,4,31 = 1 +- K— = ^ 

2+- r 

4+ - 
3 

C 5 = [1,2,4,3,2] = 1 +- — - 

2+-1— 


At first sight it appears that the only way of evaluating an FCF is the 
logical one of starting with the final partial quotient and working our way 
upward through the expression, determining the various fractions along the 
way. Fortunately this is not so. We can evaluate the convergents 
Cl, Ci, ... directly, obtaining each one from the convergents which have 
gone before it until we reach the one we want. Our next result will show 
how this is done. 


Although, for example, we have, 
written C 4 = [ 1 , 2 ,4,3] the 
convergent C 4 is really the value of 
42 

this expression, namely —. 


A key step in the proof of the coming theorem is going to be the following 
identity. 


Second continued fraction identity 


[cii, a 2,...,ctfc, o/j-i-i] — [tti, ci 2 ) • • •,o/j_ 1, o,k H-] 

Ofc+l 


This is obtained by absorbing the final partial quotient of the left-hand side 
into the penultimate one. This enables us to write a continued fraction with 
k + 1 partial quotients in an alternative form with just k partial quotients. 
However, the expression on the right-hand side of this equality is not a 

simple continued fraction because, in general, a*, -I--— is not an integer. 

Q*k +1 

For this reason we shall state, and prove, our theorem for the case of the 


a,k -I-will be an integer only 

when dfc+i = 1 . 
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general continued fraction, although thereafter our applications will be 
restricted to the special case where the continued fraction is simple. 


Theorem 1.2 The value of a finite continued fraction 

Let [oi, 02 ,..., a n ] be any finite continued fraction. Define the 
numerators Pi, P2, P3, ■■■ > Pn and the denominators 91, 92, 93, ..., q n 
of this continued fraction by: 

Pi = fli, p 2 = oio 2 + l; Pk = a,kPk-i + Pk-2, 3 <k<n; 

9 i = 1, 92 = a 25 9 fc = a k qk -1 + 9 k-2» 3 < k < n. 

Pn 

Then the value of [oi, a 2) ..., a n ] is —. 

Qn 


Before we tackle the proof, let us get a feel for the formulae involved here by 
working through some examples. 


Example 1.2 

Use the formulae of Theorem 1.2 to determine the value of [3,3,1,2,5]. 

The first two numerators are 

Pi = oi = 3 and p 2 = a i a 2 + 1 = 3x3+1 = 10. 

Subsequent numerators are given by the iterative formula 
Pk = CLkPk -1 +Pk- 2 - 

P3 = 1 x p 2 + Pi = 1 x 10 + 3 = 13 

Pi = 2 x p 3 + P 2 = 2 x 13 + 10 = 36 

P 5 = 5 x pi + P 3 = 5 x 36 + 13 = 193 

The first two denominators are 
qi = 1 and 92 = 02 = 3. 

Subsequent denominators are given by the same iterative formula as for the 
numerators, qk = afc9fc -1 + 9 fc- 2 - 

93 = 1x92 + 91 = 1x3+1 = 4 
qi = 2 x 93 + 92 = 2x4 + 3=11 
95 = 5 x 94 + 93 = 5 x 11 + 4 = 59 

The value of [3,3,1,2,5] is, according to the theorem, — = —♦ 


Notice that the calculations in the above example give us more than just the 
value of the FCF; we have determined the values of all the convergents along 
the way. For example, the theorem tells us that the third convergent [3,3,1] 

Indeed the convergents are precisely —, —, —, — and 

9i 92 93 94 

10 13 36 , 193 , , 

—, —, — and -, the last of which is the value of the 

3 ’ 4 11 59 ’ 


has value — = ~. 

4 


93 

P5 T , , . 3 

—. I hat is, 

95 1 

FCF itself. 


This example illustrates an immediate corollary to the above theorem. 


Corollary to Theorem 1.2 


Let [ 01 , 02 , ..., a„] be any finite continued fraction and the numerators 

Pi,P2, ■ ■ - Pn and denominators 91 , 92 , • • 
theorem. Then its convergents are 

., q n be as defined in the 

Ck = —, for k = 1 , 2 ,... ,n. 


9fc 
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Problem 1.5 


Use the corollary to Theorem 1.2 to determine the convergents of the 
following FCFs. 

(a) [2,3,4,1,2] (b) [-2,3,1,6] 


In the next example we suggest a tabular method for calculating the 
successive numerators and denominators. This method may help you avoid 
arithmetic slips and might well assist you in remembering these important 
formulae. However, after this subsection we shall not display tables again; 
we shall just use the formulae to write down the convergents. 


Example 1.3 

Determine the rational number represented by [4, 2, 1, 3, 1, 2, 4], 

In the table we set out four rows labelled k, pk, ak and qk . Since t his 
continued fraction has seven partial quotients, k takes the values from 
1 to 7. The ak are the given partial quotients and so the entries for this row 
can be displayed immediately. The pk and qk rows are started by entering 
the first two values for each using the formulae p\ = ai, p 2 = a\a 2 + 1, 
qi = 1 and q 2 = a 2 . 


k 

1 

2 

3 

4 

5 

6 

7 

Pk 

4 

9 






O-k 

4 

2 

1 

3 

1 

2 

4 

Qk 

1 

2 







To complete the table we move along the ak row filling in entries above and 
below each a*. For the entry above a k (i.e. pk), we multiply ak by the entry 
diagonally above and to the left, and add to this product the entry to the 
left of the latter (Figure 1.1). So, for example, p 3 = 1x9 + 4 = 13 and 
Pi = 3 x 13 + 9 = 48. 










qk -2 

- y x 





qk — akqk -1 + qk -2 


Figure 1.1 Figure 1.2 

The same procedure applies for the qk . We multiply a* by the entry 
diagonally below and to the left and add on to this product the entry to the 
left of the latter (Figure 1.2). So, q 3 = 1 x 2 + 1 = 3, and so on. 





The complete table is 


k 

1 

2 

3 

4 

5 

6 

7 

Pk 

4 

9 

13 

48 

61 

170 

741 

a-k 

4 

2 

1 

3 

1 

2 

4 

qk 

1 

2 

3 

11 

14 

39 

170 


, , 4 9 13 48 61 170 , 741 _ _ , 

and the convergents are —, —, —, —, and the hnal one being 

the value of the FCF. ♦ 

Problem 1.6 _ 

Use the tabular method to determine the convergents of each of the 
following FCFs. 

(a) [-4,1,5,1,6] (b) [1,1,2,1,2,1,2] 

Problem 1.7 _ 

In Problem 1.5 part (a) we determined the convergents of [2,3,4,1,2]. Now 
determine the convergents of the following FCFs. 

(a) [2,3,4,1,2,2] (b) [2,3,4,1,2,100] 


Problem 1.7 highlights another property of continued fractions. The 

r , 2 7 30 37 , 104 . 

convergents of 2,3,4,1,2 are -, -, —, — and ——, and any continued 

1 3 13 16 45 

fraction which begins [2,3,4,1,2,...] will have these five as its first, second, 
third, fourth and fifth convergents no matter what follows. There is nothing 
special about [2,3,4,1,2] here. If [ai, a 2 ,..., a n ] is any continued fraction 

with convergents —, —, ..., —, then the continued fraction obtained by 

<7i q-i q n 

adding a further partial quotient, [ai, a ?,..., a n , a n+ i], has the same 
convergents —, —, ..., — but with one extra final one, which will be equal 

qi <72 q n 

to the value of the (new) continued fraction itself. According to 

Theorem 1.2, this final convergent is Qn+1 ^ w ~*~^ n 1 , and this is what we 

Qn+lQn + q-n —1 

have to establish for a proof by induction. 

With this observation we are now ready to prove Theorem 1.2. 


Proof of Theorem 1.2 

We prove the result by Mathematical Induction on the number n of partial 
quotients in the continued fraction. 


As the basis for the induction, we note that the result holds for any 
continued fraction [oq] with a single partial quotient and for any continued 
fraction [ai, <22] with two partial quotients since 


M = ai 


and 


Pi 

<Zi 



and 


[01,02] = ai + — and — 
a 2 <72 


aia 2 + 1 ,1 

-— 01 H-. 

a 2 a 2 
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Now suppose that the result holds for any continued fraction with k > 2 
partial quotients and consider the finite continued fraction 

[ai,a 2 ,... ,a*;+i] 

with k + 1 partial quotients. The numerators of this continued fraction are 
Pi> P 2 , ■ ■■, Pk and Pk+i, and the denominators are q\, q 2 , .... Qk and Qk+i- 

Our task is to show that the value of the continued fraction is Pk+1 . 

Qk +1 

Now, by the second continued fraction identity, 

[«ii • • •,flfc+i] = [a 1 ,o 2 ,..., ak-i,a' k ], where a' k = ak H-—. 

flfc+i 

Look at the right-hand continued fraction. Its numerators are 
Pi,P 2 ,---,Pk-i and p k =(a k -1-— ) p fc _i + p fc _ 2 . 

V a k +\J 

Similarly its denominators are 

9i,92 ,--.,qk-i and q' k = ( a k + —) q k ^ + q k _ 2 - 

V a fe+i/ 

Moreover, as it has just k partial quotients, the induction hypothesis tells us 

that its value is . That is, 

9fc 

Pk 


[ai,a 2 , • • • ,a&+i] = — 
% 


ak + 


« fc +1 


Pk -1 + Pk -2 


( a k H-) qk -1 + qk -2 

V a k+l) 

= Qfc+i(flfcPfc-i + Pk-2) + Pk-i 
ak+i(a k q k -i + qk-2) + 9fe-i 
_ «fc+iPfc +Pfc-i _ Pfc+i 

a-k+iqk + qk -1 qk+i 

This completes the induction step and therefore, by Mathematical 
Induction, the result is true for a continued fraction with any number of 
partial quotients. 


1.3 Properties of the convergents 

The convergents of an FCF have a number of curious properties. For 
example, let us look again at the sequence of convergents of [4, 2,1,3,1,2,4] 
which we found in Example 1.3. 

Pk 4 9 13 48 61 170 74]” 

9fe 1 2 3 11 14 39 170 


The first thing you may have noticed is that, for each k, the numbers pk and 
q k are relatively prime. There is also something curious about the 


9 13 

‘cross-products’ of consecutive convergents. For example, for - and — we 
have 9 x 3 = 27 and 13 x 2 = 26. They differ by 1, something which 
happens with every pair of consecutive convergents in this example. You 
might like to check that 170 x 14 — 61 x 39 = 1. This is a general property 
of the convergents of an FCF. 





Theorem 1.3 Properties of the convergents 


Let [ai, a 2 ,..., a n ] be any FCF, and let pi, p 2 , ... , p n and 
qi, q 2 , ■ ■ ■, Qn be its numerators and denominators as defined in 
Theorem 1.2. Then we have the following properties. 

(a) PkQk-i -Pk-iQk = (~l) fc ; that is, 


Pk _ Pk -1 
Qk qk- l 


qkqk-i ’ 


for 2 < k < n. 


(b) Pkqk- 2 ~Pk- 2 qk = (-l) fc 1 ak ; that is, 


Pk _ Pk-2 
Qk Qk -2 


QkQk-2 


for 3 < k < n. 


(c) pk and qk are relatively prime for k = 1, 2, ..., n. 


Proof of Theorem 1.3 

(a) We use Mathematical Induction to prove that, for any integer r > 2 the 
numerators p r -i and p r and the denominators q r -\ and q r of the FCF 
[flj, CL 2, . . . , CLf ] satisfy 

p r q r -l - Pr-iqr = (-I) 7 "- 

For the case r = 2, p\ = ap 2 = ai a 2 + 1, q\ = 1 and </ 2 = n 2 so that 

P 2 qi ~ Pi92 = («i« 2 + 1) x 1 - c*ia 2 = 1 = (-1) 2 . 

This gives the basis for the induction. 

For the induction step, assume that it is true for r = k, so that, for the 
FCF [ai,a 2 ,...,ak}, we have that Pkqk -1 ~ Pk-iqk = and 

consider [a\, a 2 ,... ,ak, a/c+i]- Its final numerator and denominator are, 
according to Theorem 1.2, 

Pfc+i = ak+iPk + Pk -1 and qk+i = a-k+iqk + qk- 1- 
Therefore 

Pk+iqk -Pkqk+l = (ak+lPk + Pk-\)Qk ~ Pki^k+iqk + Qk- 1) 

= Pk-\qk - Pkqk -1 
= (~l)(pkQk-l ~ Pk-lQk) 

= (—1)(—l) fc , by the induction hypothesis, 

= (-l) fc+1 , 

showing property (a) to be true when r = k + 1. This completes the 
proof. 

(b) This is a simple consequence of property (a). 

PkQk-2 - Pk-2Qk = (^kPk-l + Pk-2)Qk-2 ~ Pk-2(o-kQk-l + qk- 2) 

= CLk(pk-lQk-2 - Pk- 2 qk-l) 

= Ofc( —l ) fe_1 

(c) This also follows from property (a). 

If d = gcd (pk, qk) then Pkqk-i ~ Pk-iQk, being an integer combination of 
Pk and qk, is a multiple of d. Therefore, by property (a), d divides 
(—l) fc , which implies d = 1. ■ 

Let us pause at this point to look at one application of FCFs. The result of 
Theorem 1.3 leads to a neat way of solving linear Diophantine equations. 


Theorem 1.2 told us that the 

convergents are —. Property (c) 
Qk 

confirms that these numbers are 
rationals, written in their lowest 
terms. 


The result of property (a) will 
follow by letting r range from 
2 to n 
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Example 1.4 

Find the solution of the linear Diophantine equation 13a: + 59 y = 10 in 
which x takes its smallest positive value. 

59 

We begin by finding the FCF of —. 

13 

59 = 4 x 13 + 7 
13 = 1 x 7 + 6 
7=1x 6 +l 


6 = 6x1 

So — = [4,1,1,6], and its convergents are y, -, 
6x9 + 5 59 

6x2 + 1 - 13 

guaranteed would occur), we observe that 


5 + 4 9 

--r = x an d 

1 + 1 2 

. From the final two convergents, (as Theorem 1.3 


59 x 2 - 9 x 13 = 1. 


Multiplying through by 10 and rearranging gives 
13 x (-90) + 59 x 20 = 10, 

showing that x = -90, y = 20 is one solution of 13a: + 59 y = 10. The general 
solution is therefore 


x = -90 + 59f, y = 20 - 13f, t G Z. 

For the smallest positive value that x can have in a solution take t = 2, to 
give the solution x = 28, y = — 6 . 4 


Problem 1.8 _ 

48 

By determining the convergents of —, find the general solution of the linear 
Diophantine equation 48a: + 19y = 4. 


Returning to Theorem 1.3, the alternative formulations of properties (a) and 
(b) give important information about the convergents of an FCF. 

Property (a) tells about the difference between successive convergents: 


C k - C k -1 


(- 1 )* ; 
Qkqk-i ’ 


for k = 2 ,3,..., n. 


(*) 


Now the denominators qk satisfy 


We could equally well have used 
19 

the convergents of — but taking 
48 

the larger number on the top is 
more efficient. 


Pk 

Remember that Ck = —■ 
Qk 


qi < qi < <73 < • • • < q n - 

Therefore the right-hand side of equation (*) becomes smaller as k increases. 
That is, the convergents are coming closer together. 

Property (b), on the other hand, tells about the difference between 
convergents two apart: 

Ck~Ck~ 2 = -— --, for k = 3,4,... ,n. (**) 

qkqk- 2 

We shall use both these facts in the proof of the following theorem 
concerning the relative sizes of convergents. For ease of notation we refer to 
the convergents C2, C4, Cq, ..., which have even subscripts, as the even 
convergents whilst Ci, C3, C5, ... are called the odd convergents. 


We ask you to verify this in 
Problem 1.10 
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Theorem 1.4 The relative sizes of convergents 

For any FCF the odd convergents form a strictly increasing sequence 
Ci < Cz < C 5 < • • •, while the even convergents form a strictly 
decreasing sequence C 2 > C 4 > Cq > ■ ■ ■. 

Every even convergent is greater than every odd convergent. 


We can illustrate pictorially what the theorem is claiming by imagining a 
plot of the points C r against r. 



• C„ 


Figure 1.3 The convergents of an FCF 

The final convergent C n , which is the value of the FCF, may be either an 
even or an odd convergent. The last part of the theorem tells us that 
whichever it is, C n is at least as large as every odd convergent and no larger 
than any even convergent. 

Proof of Theorem 1.4 

As > 0 (for k > 2) and qk> 0 (for k > 1) we note that, for k > 3, 

QkQk—2 

has the same sign as (—l ) fc_1 and so will be positive when k is odd and 
negative when k is even. 

Putting k = 2r + 1 in equation (**) produces 

C 2 r+i - C 2r -i > 0 , for r > 1 , 

which confirms that the odd convergents are strictly increasing. 

Putting k = 2r in the same equation gives 
C 2r - C 2r - 2 < 0 , for r > 1, 

which confirms that the even convergents are strictly decreasing. 
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It remains to show that every odd convergent is less than every even 
convergent, and for that we use equation (*) which tells us that — Ck -1 

takes the same sign as (-l) fc , for k = 2, 3 , ..., n. Putting k = 2r gives 

C 2r — C2r-1 > 0, for T > 1, 

and putting k = 2r + 1 gives 

C*2r+1 — Cir < 0, for r > 1. 

Hence the even convergent C 2r is larger than each of the adjacent odd 
convergents C 2r -i and C2 r +i- 

Now consider any even convergent C 2s and any odd convergent C 2t+ i • 

If s >t then C 2s > C 2 s+i > C 2 t +i, the second inequality utilizing the fact 
that the odd convergents increase. 

If s < t then C 2s > C 2t > C 2 t + 1, the first inequality utilizing the fact that 
the even convergents decrease. 

In each case C 2s > C 2t +\, as required. ■ 

Problem 1.9 _ 

Determine the convergents of [2, 1, 1, 2, 1, 1, 2] and confirm that they 
satisfy Theorem 1.4. 

Problem 1.10 _ 

Show that, for any FCF, the denominators qk are increasing, and are strictly 
increasing with the possible exception of qi = q 2 . Show further that, for any 
FCF, qk> k for all k > 5. 





2 INFINITE CONTINUED FRACTIONS 


2.1 The convergence of a continued fraction 

In the introduction to the unit we built up a sequence of FCFs 

[ 1 ], [ 1 , 1 ], [ 1 , 1 , 1 ], [ 1 , 1 , 1 , 1 ],... 

whose rational values give better and better approximations to the golden 

ratio, However, the value of any FCF is rational, but the golden 

ratio is irrational, and so no matter how far we go along this sequence of 
FCFs we can never reach the golden ratio itself. What is of interest is the 
limit of this sequence of rationals. This suggests that we might explore the 
idea of the unending continued fraction [1,1,1,...]. 

This is our first example of a simple infinite continued fraction. 


Definition 2.1 Infinite Continued Fractions. 

Let a\ be any real number and let «2, 03, ..., a n , ... , be any positive 
real numbers. The expression 

1 

01 H- 1 - 

h2 H- 1 - 

a3 H- 

'■ r 

a n -\ - 

is an infinite continued fraction. We write this continued fraction in 
terms of its partial quotients as [a,\ , 0,2,0.3,.. .]. Its convergents are the 
finite continued fractions 

Ci = [ffli], C2 = [01,02], C3 = [01,02,03],.... 

When all the partial quotients are integers, the continued fraction is 
said to be simple. We shall use the abbreviation ICF for a simple 
infinite continued fraction. 


In this subsection we shall explore basic properties of infinite continued 
fractions. Our ultimate interest will lie with simple infinite continued 
fractions although, as in the previous section, we shall have occasion to allow 
the partial quotients to be non-integers for the purpose of providing general 
proofs. 

What exactly does an infinite continued fraction represent? Whereas 
calculating the value of a finite continued fraction involves only the normal 
operations of arithmetic, evaluation of an ICF is an infinite process and as 
such requires a different approach. There is a strong resemblance to the 
convergence question for infinite series. The sum of the infinite series 

OO 

^2 X r = X\ + X2 + X 3 -)- 

r= 1 

is defined to be the limit of the sequence 

an, an+X 2 , X 1 +X 2 +X 3 , ..., X1 +X 2 +X 3-I- \-x n , ... 

of partial sums, provided this limit exists. For ICFs the convergents fulfil 
the role of the partial sums and so, by analogy, we define the value of an 
ICF to be the limit (if it exists) of the sequence {C„} of its convergents. 


We use the notation {C„} to 
represent the infinite sequence 

Ci,C 2 ,C 3 ,... . 
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Definition 2.2 Value of an Infinite Continued Fraction 

The value of the infinite continued fraction [ai, a 2 , a 3 ,...] is defined to 
be the limit of its sequence of convergents {C„}. 


Notice that we did not add the condition ‘provided this limit exists’. This is 
because there is one big difference between infinite series and ICFs: whereas 
some infinite series converge and some do not, it transpires that for every 
ICF, its sequence of convergents has a limit. 


Theorem 2.1 Sequences of convergents converge 

For any ICF [c*i, a?, 03 ,...] its sequence of convergents {C„} tends to a 
limit. 


Proof of Theorem 2.1 

For any n > 1, the the first n convergents of the ICF [ai, a 2 , a 3 ,...] are 
Ci, C2, ..., C n - 1, C n . These are precisely the convergents of the FCF 

Pfl 

Cn = [«i,«2, 0.3,..., a n ]. What is more, C n = —, where p n and q n are the 

Qn 

numerators and denominators as determined by the formulae in 
Theorem 1.2. 

Now Theorem 1.4 tells us that 

Cl <C 3 <--< C 2 k+\ < ■•• <C n < •■• < C 2 k < • • • < C4 < C2, 

where 2k + 1 < n. As such a chain of inequalities holds for any value of n, 
no matter how large, we see that the odd convergents of the ICF form a 
sequence Ci, C 3 , C 5 , ... which is strictly increasing and bounded above 
by C2. A well-known theorem of analysis tells us that a strictly increasing 
sequence of real numbers which is bounded above is convergent, so 
Ci, C 3 , C5, ... has a limit, say a. 

In the same way the even convergents form a strictly decreasing sequence 
C2, C4, Ce, ... which is bounded below by Ci, and so it has a limit, say b. 


\ 


nC 2 






Figure 2.1 



We shall use the term ‘converges’ 
as an alternative to ‘tends to a 
limit’. So the following theorem 
tells us that the sequence of 
convergents of every ICF converges! 


This result usually goes by the 
name of The Monotone 
Convergence Theorem. Here, we 
are going to assume this result is 
known, along with some other 
results on the convergence of 
sequences of real numbers. 
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For the whole sequence C\, C 2 , C3, C4, ... to converge it remains to show 
that a = b. As every even convergent is greater than every odd convergent, a 
and b both lie between C^n-i and C 2n , for each positive integer n. Now 


C 2n ~ C 2n -i = -> by property (a) of Theorem 1.3, 

Q2nQ2n—l 


_ 1 _ 

(2n)(2n — 1) ’ 


since qk > k for large k. 


Therefore 

|a - 6| < C 2n ~ C 2n -1 < 77 TT77, - TTi 

(2n)(2n — 1) 

and, as the term on the right tends to 0 as n increases without bound, we 
conclude that a = b. ■ 


Problem 1.10. 


The above proof does more than confirm that the limit exists; it contains an 
estimate for the accuracy of the convergent C n as an approximation to the 
actual limit. We know that the limit lies somewhere between any pair of 
consecutive convergents. If we couple this with the information given by 
Theorem 1.3 property (a), which gives a bound for how far apart successive 
convergents can be, we have the following immediate consequence. 


Corollary to Theorem 2.1 

If {C n } is the sequence of convergents of x then, for each n > 1: 
(a) x lies between C n and C' n+ 1; 

1 

< -. 

QnQn+l 


(b) 


Pn 

Qn 


When we talk of the convergents 
of x we mean the convergents of 
the ICF of x. 


Pn 

< 

Pn+l 

Pn 

1 

x — — 




Qn 


Qn+1 

Qn 

QnQn -\-1 


Problem 2.1 - 

Estimate the maximum possible error in taking the convergent C5 as an 
approximation to the value of [1,1,1,1,...]. What is the maximum error in 
C 6 ? The actual value of this ICF is the golden ratio. Taking the latter to be 
1.61803 to five decimal places, what are the actual errors in C5 and C $? 


2.2 Periodic continued fractions 

We now know that any ICF has a value, but how do we go about finding that 
value? One thing we can do is determine more and more of its convergents 
in the knowledge that we are obtaining better and better rational 
approximations to its value. More about that later, but for the time being 
there are some ICFs whose value we can find without too much difficulty. 


Example 2.1 

Find the value of the ICF [1,3,2,1,3,2,1,3,2,...], where the sequence 1,3,2 
of partial quotients is repeated indefinitely. 


Knowing that the ICF converges, let its value be x. Then notice that 


1 + 


1 


= 1 + 


1 


1 + 


3 + 


3 + 


3 + 


1 


= [l,3,2,x], 


2 + 


2 + 


1 + 


[1,3,2,...] 


2 + — 
x 


3 + 


2 + — 
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and, as this continued fraction has value x, 
x = [1,3, 2, a;]. 

The convergents of [1,3,2, x] are, using the formulae for the numerators and 
denominators, 

1 4 9 9a; + 4 

1’ 3’ 7 7x + 3 

and, as the final convergent is x itself, we have the equation 

9a; + 4 9 

x = --that is 7ar — 6 a; — 4 = 0. 

7x + 3 

Clearly x > a\ = 1 and so we want the positive root of this equation. Using 
the formula for the roots of a quadratic equation we obtain x = ---. ♦ 

Problem 2.2 _ 

Consider the ICF [3,1,3,3,1,3,3,1,3,...], where the sequence of partial 
quotients 3, 1, 3, recurs indefinitely. Determine its value x by writing 

x = 3 H-j— 

1 + ” r 

3 H— 
x 

and solving the resulting quadratic equation. 

Problem 2.3 _ 

Find the value of the ICF [2,2,2,2,2,...]. From this result write down the 

ICF for \[2 and hence find a rational number — which lies within —-— 

b 1000 

of y/2. 

In Example 2.1, where we effectively wrote 
x = 1 +-—, 


we were taking some liberties with properties of limits. To give the required 
argument correctly we should really be dealing with finite expressions, and 
subsequently taking limits in the knowledge that they do exist. That is, 
instead of working with the ICF whose value is x we should deal with its nth 
convergent, C n . From our work on FCFs we can certainly write 


C n — 1 + 


for n > 4, 


which leads (as above) to 

9C»-3 + 4 
” 7C n _ 3 + 3 

Now as n increases without bound we know that C n converges to x, and 
consequently C n _3 converges to x. Taking limits of both sides gives 

9x + 4 

x =-, 

7x + 3 

which is the situation we reached in Example 2 . 1 . 
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The property which we have been using is the following. 


Infinite continued fraction identity 

— 


r i 1 

[ai, a 2 ,d3,.. .J — dj + .. 

[d 2 ,d3,d4,...J 




This follows since we have the corresponding result for FCFs 


[fllj 0,2, ■ ■ ■ , ®n] — Ol + 



and, as n increases without bound, the FCFs on each side of this equation 
converge to the required ICFs. 

So far, the only ICFs which we have been able to evaluate are ones in which 
a block of partial quotients repeat indefinitely. We have a name and 
notation for such ICFs. 


Definition 2.3 Periodic Continued Fractions 

An ICF which consists of a block of partial quotients which repeat 
indefinitely is said to be purely periodic, and we write it as 
[(ai, d 2 ,0.3 ,..., a n )] enclosing the repeating block in angular brackets. 
We refer to the repeating block as the cycle of the ICF. 

An ICF of the form [fq, b 2 ,..., b m , (di, d 2 , ■ ■ ■, a n )], in which a 
repeating block occurs after a finite number of partial quotients, is said 
to be periodic. 


For example, 

[3,4,3,4,3,4,...] = [(3,4)] 

is a purely periodic ICF, and 

[1,1,3,4,3,4,3,4,3,4,...] = [1,1, <3,4)] 

is periodic. Periodic continued fractions may be written in different forms 
using the angle bracket notation. For example, 

[(3,4)] = [3,(4,3)] = [3,4, (3,4)]. 

The evaluations that we have carried out illustrate a general strategy for 
finding the value of a purely periodic continued fraction. To summarize, 
consider [(di, d 2 , <13, ..., a n )\ with value x. The key step is the observation 
that x can be written as the finite (but not simple) continued fraction 

X — [di , 0 .2 , 0-3) • • • > 1 x\ • 

Adopting the usual notation for the convergents, the final two convergents of 
Pn 1 Pfi 

[dj, a 2 , 03, ..., a„l are yn ~ and —, and can be determined by the methods 
q n - 1 Qn 

of the previous section. The final convergent of [di, a 2 , 03, ..., a„, x] is then 

XJpYi + Pn _ 1 

— -—, and equating this to x yields the quadratic equation 

xq n + q n -1 

q n x 2 + (q n - 1 -Pn)x-p n - 1 = 0. 

As the first partial quotient ai reappears infinitely often as a partial 
quotient it has to satisfy a\ > 1, and so the root we seek is certainly 
positive. In fact this quadratic has one positive and one negative root since 

the product of its roots is equal to-which is negative since p n ~\ and If ax 2 + bx + c = 0 has roots a and 

q n are each positive. P then a0 = - and a + /3 = — 
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It is only a short step from the evaluation of purely periodic continued 
fractions to the evaluation of any periodic continued fraction, as illustrated 
by the following example. 


Example 2.2 

Evaluate the periodic continued fraction [—2,1, (3,4)]. 

Let y = [-2,1, (3,4)] and x = [(3,4)], so that y = [-2,1, x]. 

We evaluate the purely periodic x by noting that x = [3,4, xl, the latter ICF 
, . 3 13 13x + 3 

having convergents -, — and ^ . Equating the final convergent to x 

gives the quadratic 

4x 2 - 12x - 3 = 0, 


and this has positive root x = 


3 + 2>/3 
2 


_2 j ^ a 

Turn now to y = [-2, l,x]. It has convergents —, — and —-, and so 

1 1 x +1 

3 + 2^3 „ 

= -x-2 _ 2 2 _ -7-2y/3 

* + 1 3 4- 2\/3 | x _ 5 + 2V3 

_ -7-2V3 5-2v/3 

5 + 2x/3 5-2v/3 

_ —23 + 4\/3 

13 ^ 

The method illustrated in Example 2.2 can be used to evaluate any periodic 
ICF. 

Problem 2.4 .__ 

Evaluate the following ICFs. 

(a) [2, (1,1,1,4)] (b) [2,2, (3)] 


Note carefully the final step in the 
calculation in this example. In 

, .. a + bJn 

general the expression - 

c + d\fn 

may be simplified by multiplying 
top and bottom by c — d\fn and 
appealing to the identity 

(c + dy/n)(c — d\/n) = <? — nd 2 . 


The net result is that all 
square-roots will be removed from 
the resulting denominator. 
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3 THE CONTINUED FRACTION 
ALGORITHM 


3.1 The continued fraction of an irrational 

Having seen that every finite simple continued fraction has a rational value 
and, conversely, that every rational number can be represented by a finite 
simple continued fraction, the following result should come as no surprise. 


Theorem 3.1 Irrationality of ICFs 

The value of any infinite simple continued fraction is an irrational 
number. 

You may have thought that this theorem must follow immediately, but we 
must exclude the possibility that a rational number could be represented by 
both a finite and an infinite simple continued fraction. 


Proof of Theorem 3.1 

Let x = [ax, a 2 , 03 ,...] be an ICF. For each n > 1, x lies between the 

consecutive convergents C n = — and C n+ i = ^ >n+ - , and so 

Qn Qn+1 


X — 


Pn 

Qn 


< 


1 


QnQn+1 


Now suppose that x is rational; say x — -, where r and s are integers with 
s > 0. Then 


0 < 


r _ Pn 
S q n 


I rq n ~ sp n \ 


< 


1 


sq n 


QnQn+1 


so that 

o < I rq n - spnl <-• 

qn+i 

But the denominators q n form a strictly increasing sequence of integers and 
so we can choose n so that q n +1 > s. But then \rq n — sp n | would be an 
integer lying strictly between 0 and 1, which is impossible. The conclusion 
from this is that no such integers r and s can exist; in other words, x must 
be irrational. • 


See the Corollary to Theorem 2.1. 


Notice that 
r ^ Pn 
S Qn 

for otherwise the continued fraction 
would be finite. 


We turn now to the converse problem. Given an irrational number can it be 
expressed as a simple infinite continued fraction? In answering this question 
in the affirmative we shall discover more: the simple continued fraction of an 
irrational number is unique. Our approach will be to demonstrate a method, 
called the Continued Fraction Algorithm, for determining an ICF for any 
given irrational number. Of course we cannot determine all its partial 
quotients unless there is some recurring pattern such as occurs with periodic 
ICFs. What the algorithm does is enable us to determine as many partial 
quotients as we please. The situation is analogous to the problem of finding 
the decimal expansion of 7r; the full expansion of tt cannot be given, but it 
can be determined to any required number of decimal places. 

Given an irrational number x\ we seek integers ai, a 2 , ... so that 

X\ [fli, fl 2 , . . .]. 

Knowing that x x lies strictly between any successive pair of its convergents, 

it has to lie between C x = a, and C 2 = a x H-. In particular, as a 2 > 1, x x 

a 2 

lies strictly between the integers a x and a x + 1. It follows that a\ = int^i), 


o 


< 
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Remember int(ai) denotes the 
integer part of x. 





and we note, in particular, that a\ is determined uniquely from X\ alone. 
Now let us write 


X\ = int(xj) + frac(xi), 

where frac(xi) denotes the fractional part of x\. Then 

1 ... 1 


X\ — di + 


a 2 


1 


int(xi) + 


a 2 + 


1 


= int(xi) + 


a 3 + 


1 

a 3 H- 


[02,03,04,. 


and so 


[02,03,04,...] 


1 

xi — int(xi) 


1 

frac(xi) 


If we write X 2 


then X 2 = [a 2 , 03 , 04 ,...] and X 2 > 1 , as it is an 


frac(xi) 

irrational number greater than the positive integer a 2 . Repetition of the 
above argument now leads to a 2 = int(x 2 ), and so 02 is determined uniquely 


from x 2 and therefore from xi. Continuing in this way, writing x 3 =- 

frac(x 2 ) 

we have x 3 = [ 03 , 04 , 05 ,...] and 03 = int(x 3 ), etc. This is the Continued 
Fraction Algorithm which is summed up in the flowchart of Figure 3 . 1 . 


The fractional part of an irrational 
number lies strictly between 0 and 
1 , and so its reciprocal is greater 
than 1. 



Figure 3.1 The Continued Fraction Algorithm 
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3.2 Using the Continued Fraction Algorithm 

Before pursuing the theory let us see the algorithm in action. 


Example 3.1 

Determine the continued fraction of y/6. 

Suppose that x\ = \/6 = [a \,(i 2 , a 3, ...] and we are to use the algorithm to 
determine these partial quotients. Since 2 2 < 6 < 3 2 , y/6 lies between 
2 and 3. Therefore int(-\/6) = 2 and y/6 = 2 + (y/6 — 2) expresses it as 

integer part, 2, plus fractional part (y/6 — 2). So ai = 2 and X 2 

where X 2 = [( 12 , 03 , • • •]• 

To find a 2 we require the integer part of x?- 


VS-2’ 


X2 


y/6 + 2 V6 + 2 


V6-2 V6-2 V6 + 2 


Now, as y/6 lies between 2 and 3, y/6 + 2 lies between 4 and 5, and so the 

_j_ 2 

integer part of —-— is 2. Therefore 


x 2 = 2 + 


2 

V6-2 


(*) 


expresses X 2 as its integer part plus its fractional part. So 02 = 2 and 
2 

x 3 = —p -• 

V6- 2 

We continue by finding the integer and fractional parts of X 3 : 

2 2 y/6 + 2 


X 3 = 


y/6-2 
So 03 = 4 and X 4 = 


y /6 - 
1 


= V6 + 2 = 4 + (V6-2). 


2 y/6 + 2 

-. But now X 4 = X 2 and we are in a cycle. We 


(**) 


y/6 — 2 

repeat calculation (*), establishing 0,4 = 2 , then repeat (**), establishing 
a 5 = 4 , then back to (*) and so on. The sequence of partial quotients 
repeats the cycle 2, 4 indefinitely from here onwards and we conclude that 
V6=[2,<2,4)]. ♦ 


Note the trick of multiplying top 
and bottom by y/6 + 2 to remove 
square roots from the denominator. 


The integer part of y/6 + 2 is 4. 


Our second example is similar, but this time we shall give just the 
calculations involved without the commentary. 


Example 3.2 


Determine the continued fraction of 


7+ y/2 


X\ 

X2 = 

X 3 = 


7+ y/2 
2 
2 


4 + 


y/2- 1 


fl! =4, X2 = 


V 2 + 1 


y /2 — 1 v /2 + 1 

1 2 V 2 + 2 _ y/2 + 1 

2 V 2-2 X 2 V 2 + 2 ~ 2 


2y/2 + 2 = 4 + ( 2 V 2 - 2); a 2 =4, x 3 = 

y/2 - 1 


1 + 


a 3 = 1, X 4 = 


y/2 — 1 
1 

2 >/ 2-2 

2 

y/2 — 1 


x'4 = X2, and so 04 = 0,2 , 0,5 = 03, — 0,4, etc. 


Therefore 

7 + y/2 
2 


[4,(4,!}]. 


♦ 
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Problem 3.1 _ 

Use the Continued Fraction Algorithm to determine ICFs for the following 
irrational numbers. 

(a) (b) 2 + V7 


Having seen illustrations of the Continued Fraction Algorithm in use, we 
must now do a little tidying up. We have seen that the algorithm will always 
produce an ICF, but can we be certain that the value of this ICF is equal to 
the irrational number from which we started? There is also a question of 
uniqueness to be resolved; is the ICF of an irrational number unique? 


Theorem 3.2 Uniqueness of the ICF 

Given any irrational number x, the ICF determined from x by the 
Continued Fraction Algorithm converges to x, and no other ICF 
converges to x. 


Proof of Theorem 3.2 

From the irrational number x the Continued Fraction Algorithm produces 
a sequence of (non-simple) finite continued fractions as follows: 

X — X\ — [fll , X 2 ] — [di , 0.2 , 3 : 3 ] — [dl 7 d2, d3, *£4] * [dl 7 d2 7 d3, • • ■ 7 d n , X71-1-1], 

where, for all n > 1, a n = int(x n ) and the step from one stage to the next is 
accomplished by 

1 

X n — d n “b * 

X'li 4 1 

To show that the ICF which is being built up in this way has limit x we 
must show that the sequence of finite continued fractions 

M, > [dr,d2,• ■ * 7 d n ], 

converges to x. 

Using the standard notation for the numerators and denominators, 

Pk 

\ai, 02, ■ ■ ■, dfcl = — 7 w e have for n > 2 
Qk 

1 

< - . 

QnQn+1 

As the denominator sequence {q n } is, for n > 2 , a strictly increasing 

1 Pn 

sequence of positive integers, -tends to 0 , confirming that — can be 

QnQn+1 Qn 

made as close as we please to x by choosing n sufficiently large. Hence the 
ICF determined from x by the algorithm does indeed converge to x. 

It remains to show that the ICF of an irrational number is unique. To that 
end suppose that the ICFs [cti, 02 , 03 , • • •] and [ 61 , 62 , 63 ,...] have the same 
value. We shall use the Second Principle of Mathematical Induction to show 
that a n = b n for all n > 1, and so the two ICFs are identical. 

First we observe that if [ai, 02 , < 23 , • • •] = [ 61 , 62 , 63 ,...] then ai = b\ since 
both are equal to the integer part of the common value. This gives the basis 
for the induction. 

For the induction step suppose that, for some k > 1, 
ai = b\, 02 —b 2 , ■ak = bk- 
Our goal is to deduce that dfc+i = b^+r- Now, 

[dl, d27 a 3, ■ ■ •] = [dl, d 2 , . . . ,dk, [ok+l,Ok+2i • • •] ] 


x — 


Pn 

Qn 


By the Corollary to Theorem 2.1. 
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and similarly 

[h,b 2 , b 3 ,...] = [ 61 , 6 2 , ■■■ ,bk, [bk+i,bk+ 2 , •••]]• 

As these two have the same value, and have the same first k partial 
quotients, it follows that 

[<*fc+i> a fc+2, • ■ •] = [6fc+i, bk+ 2 , ■ • •]• 

But now dk+i = bk +1 as each is equal to the integer part of the value of this 
ICF, and the induction step is complete. Hence the two ICFs having the 
same value are the same. ■ 


3.3 The continued fraction of y/n 

Let us take stock of where we have reached with simple continued fractions. 
For FCFs, the value of any FCF is rational and conversely every rational can 
be represented as an FCF. However, the representation is not unique; indeed, 
each rational number can be represented as an FCF in exactly two ways. 

Turning to ICFs, the value of any ICF is irrational, and conversely every 
irrational number can be represented as a ICF. But this time the 
representation is unique; every irrational number has a unique ICF. 

The only ICFs which we have been able to evaluate so far have been 
periodic ones, and these have always turned out to involve square roots in 
some way. We are not going to digress to prove it here, but it turns out that 
the values of periodic continued fractions are precisely the numbers of the 
form r + Sy/n, where r and s are rational and n is a non-square positive 
integer. That is, each periodic ICF has value equal to one of these numbers 
and, conversely, each such number has a periodic ICF. For the present, let us 
investigate further the ICF of certain square roots. 

Example 3.3 

For each positive integer n determine the ICF for y/n 2 + 1. 

The integer part of y/n 2 + 1 is n, and so we write 
y/n 2 + 1 = n + ^y/n 2 + 1 — n'j . 

This gives 0,1 = n. We look next at the reciprocal of y/n 2 + 1 — n, which is 
simplified by multiplying top and bottom by y/n 2 + 1 + n. 

— ; -= y/n 2 + 1 + n = 2n+ (y/n 2 + 1 — n) . 

y/n 2 + 1 — n V / 

So 0.2 = 2 n, and the same fractional part y/n 2 + 1 — n has arisen. This 
means that we are immediately in a cycle and so y/n 2 + 1 = [n, (2n)]. 

For example, for n = 1, 2, 3 and 4 we obtain, respectively, the following. 
v/2 = [l,<2>] 
v / 5 = [2,<4)] 
v'To = [3, (6)] 

VT7 = [4, (8)] ♦ 


-This proof is formalizing the 
repeated use of the infinite 
continued fraction identity. From 

[<* 1 , 02 , 03 , ... ] = [61,62, f> 3 ) • • ■ ] we 

deduce that ai = b\ and 
[02,03,(14,...] = [62,1*3,64, • ■ •]• 
But then a2 = 62 and 
[03,04,...] = [63,64,...] etc. 


These numbers are called the 
quadratic irrationals and arise as 
roots of quadratic equations 
ax 2 + bx + c = 0, where o, b and c 
are integers. 
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Problem 3.2 ___ 

Determine the continued fraction of y/n 2 + 2 for any integer n> 1. 

Hence find ICFs for, y/3, y/6, n/TT and -\/l8. 


By now you should be feeling confident that you could use the Continued 
Fraction Algorithm to determine the ICF of y/n, for any positive integer n. 
In the table below we collate some of the ICFs of y/n which we have seen 
and give several more. 

Table 3.1 Simple continued fractions for y/n 


y/2 — [1, <2>] 

^3 =[1,(1,2)] 
v/5 =[2,(4>] 

y/& =[2,(2,4)] 

V7 =[ 2 ,( 1 , 1 , 1 , 4 )] 
\/8 =[2,(1,4)] 

\/IO = [3, (6)] 

\/lT = [3, (3,6)] 

Vl2 = [3, (2,6)] 


yi3 = [3, (1,1,1,1,6)] 
^14 = [3, (1,2,1,6)] 

VT5 = [3, (1,6)] 
VT7=[4,(8)] 

\/l8 = [4, (4,8)] 
yi9 = [4, (2,1,3,1,2,8)] 
v/20=[4, (2,8)] 
V2l=[4, (1,1,2,1,1,8)] 
y/n = [4,(1,2,4,2,1,8)] 


v/23 = [4,(l,3,l,8)] 

'M = [4, (1,8)] 

\/26 = [5, (10)] 

\/27 = [5, (5,10)] 

\/28 = [5, (3,2,3,10)] 

V29 = [5, (2,1,1,2,10)] 

V30 = [5, (2,10)] 

n/31 = [5, (1,1,3,5,3,1,1,10)] 

V32 = [5,(1,1,1,10)] 


y/46 =[ 6 ,( 1 , 3 , 1 , 1 , 2 , 6 , 2 , 1 , 1 , 3 , 1 , 12 )] 

y/61 = [7, (1,4,3,1,2,2,1,3,4,1,14)] 

y/M = [9, (1,2,3,1,1,5,1,8,1,5,1,1,3,2,1,18)] 

V97 =[9,(1,5,1,1,1,1,1,1,5,1,18)] 

y/l51 = [12, (3,2,7,1,3,4,1,1,1,11,1,1,1,4,3,1,7,2,3,24)] 

%/l66 = [12, (1,7,1,1,1,2,4,1,3,2,12,2,3,1,4,2,1,1,1,7,1,24)] 

Problem 3.3 ___ 

Examine Table 3.1. 

(a) What common features do you observe in the continued fractions? 

(b) Can you draw any conclusions concerning the continued fraction of 

y/n + int(v^)? 

(c) Conjecture the general form of the continued fraction of y/n, for n a 
positive non-square integer. 


From our solution to Problem 3.3 we are led to believe that, for any 
non-square positive integer n, the ICF of y/n takes the form 

y/n = [ai, ( 02 , 03 ,... , a 3,fl2,2oi)], 

where a\ is the integer part of y/n. Here the periodic part may consist of: 

• just (2ai) as in, for example, y/2 = [1, (2)] and y/l7 = [4, (8)]; 

• an odd number of partial quotients as in y/l3 = [3, (1,1,1,1,6)]; 

• an even number of partial quotients as in y/l = [2, (1,1,1,4)]. 

It transpires that the general form is, indeed, as claimed. However proof of 
this would take us deep into the study of more properties of continued 
fractions, which we do not have space for in this course. 




4 RATIONAL APPROXIMATIONS 


4.1 Approximating irrationals by rationals 

We have seen that an irrational number has a unique ICF and that the 
convergents provide a sequence of better and better rational approximations 
to the irrational number. We even have a way of estimating how good these 
rational approximations are. In this section we shall go further and show 
that continued fractions provide what are, in a certain sense, best rational 
approximations. 


Given an irrational number x, we can always improve on any rational 
approximation - to x. In fact there are infinitely many rational numbers 

lying between - and x. Consequently there is no such thing as the best 
rational approximation to x. However, if we impose restrictions on the size 
of the denominator we can ask about best approximations. For example, of 
all rational numbers whose (positive) denominator does not exceed 50, the 

22 

best rational approximation to n is —. We have to allow the denominator 

179 

to climb to 57 before we find the better approximation-to n. 

57 

X) 

What we shall show in this section is that if — is a convergent of x then 

CL ^ n 

there is no rational number — with 0 < b < q n which approximates x better 
Pn 

than —, and in this sense the convergents of x are best rational 
qn 

approximations to it. But before becoming involved in the analysis we give 
an elegant geometric interpretation of the convergents to a continued fraction 
of an irrational number which was first given by Klein in 1895. It illustrates 
precisely the sense in which convergents are best rational approximations. 


Remember that when we write 

for a rational number we assume 
that a, b are integers with b 
positive and gcd (a,b) = 1. 


Let a > 0 be irrational and consider the part of the line y = ax in the first 
quadrant of the Cartesian plane as illustrated in Figure 4.1. This line cannot 
pass through any of the lattice points (i.e points whose coordinates are both 
integers) except (0,0), for if it did a would be rational. Imagine that p ins 
are set upright at each of the lattice points and a fine string is placed along 
the line with one end fixed to the line at an infinitely remote point, and the 
other end at the origin. Keeping the string tight we move the end at the 
origin towards the right. It will catch on a series of pins below the line. In 
fact, the pins on which it catches are those at lattice points with coordinates 
(<7i>Pi)> (q:i,P3), ( 95 ,^ 5 ), • • • corresponding to the odd convergents 
Pl_ P3 P5 

91 ’ 93 ’ 9s ’ ” ’ 

of a. Similarly, if we move the string from the origin to the left it catches on 
pegs at (r/ 2 ,p 2 ), (94>P4), (%,Pe), • • •, corresponding to the even convergents 
P2 Pi Pe 

92 ’ 94 ’ 96 ’ 

of a. 


If y = ax passes through (m, n ) 

then a = —. 

m 


Recall that 

Ci < C 3 < C 5 < ■ ■ ■ < a < • • • < C 4 < C 2 , 

and notice that Ci = — give the 
<7i 

gradients of the lines joining the 
origin (0,0) to the points ( qi,Pi ). 
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Figure 4-1 Convergents as best approximations 
to the golden ratio 


Figure 4.1 illustrates all this for the case a = the golden ratio, 

, 1 2 3 5 8 13 ^ , , , , . , 

whose convergents begin . ihe dashed line is the 

1 1 2 3 5 8 

portion of y = ax in the first quadrant. Moving to the right it catches on 
pins with coordinates (1,1), (2,3), (5,8),..., whilst moving left it catches on 
the pins at (1,2), (3,5), (8,13), and so on. 


Let us return to the question of best rational approximations. First we must 
prove, formally, a result which we have seen illustrated in the many 
sequences of convergents that we have calculated; each convergent is a better 
approximation to the irrational number than is its predecessor. The theorem 
which follows includes this result, but also gives some bounds for the 
accuracy of the convergents. 
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Proof of Theorem 4.1 

Let x have ICF [ai, < 12 , 03 ,...]. The Continued Fraction Algorithm yields the 
system of equations 

X [^15 X2] [®1 > ^2 > X3] ‘ 1 0,2 j • • • j x„ +1 ] — • * • 

and 


Pn 

Qn 


[01, • • • j *^n-|-l] 

Qn 


Xn+lPn H“ Pn— 1 


Xn+lQn Qn— 1 
Pn—lQn ~ PnQn —1 


Pn 

Qn 


by Theorem 1.2, 


9n(^n+l9n H“ Qn— l) 

= —7 ---r, by Theorem 1.3 property (a). 

Qn(x n +iq n + Qn-l) JV ' 

NOW X n _j_i = a n +2, a n-t-3? • • •] SO < X n _|_i < H-1. 


Therefore 


Pn ^ 1 1 

£ — - _ = _ 

Qn Qni.&n+lQn T Qn— l) QnQn+1 
which gives the left-hand inequality when n = k + 1 . But we also have 


_1_ 

Qn((o-n +1 + 1 )Qn + Qn-l) 

_1_ _ 1 

Qn{ a n+lQn + Qn-l) + Qn Qn{Qn+l + Qn) 
1 _ 1 

Qn{Qn+ 1 T<?n+l) 2 QnQn+1 

which gives the right-hand inequality by taking n = k. 


x — 


Pn 

Qn 


Remember that {q n \, for n > 2, is 
a strictly increasing sequence. 


For the middle inequality note that 


Qk+ 2 — dfc+29fe+i + Qk > Qk + Qk = 2 qk 
and so 


< 


Qk+iQk+2 2 qk+iqk 


The left-hand inequality in Theorem 4.1 gives an indication of how close 
each convergent is to the true value whilst the right-hand inequality gives a 
bound which limits the accuracy. We can collate the two facts for each 
convergent as follows. 


Corollary to Theorem 4.1 


Pk 


For all k > 1 the convergent — of x satisfies 

Qk 


2<Mfc+i 


< 


x — 


Pk 

Qk 


< 


1 


Qkqk+i 


Problem 4.1 __ 

The ICF of \/5 is [2, (4)]. Estimate the accuracy of the third convergent as 
an approximation to \/b. 
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4.2 Best approximations 

The discovery that the ICFs of sjn can be described exactly in terms of 
finite recurring cycles is an interesting fact, because the decimal 
representations of these numbers (when n is non-square) have no discernible 
pattern. Of course ICFs may have other patterns besides recurring blocks 
which permit them to be described exactly. Amazingly, one of the 
best-known of all irrational numbers, e, the base of the natural logarithms, 
has such a pattern to its ICF. Euler discovered in 1737 that 

e = [2,1,2,1,1,4,1,1,6,1,1,8,1,1,10,1,...], 

where the even positive integers 2, 4, 6, ... occur in order as every third 
partial quotient, the other partial quotients all being 1 with the exception of 

tti — 2. 


Example 4.1 

Determine the first convergent of e which is accurate to within 0.0001 of the 
true value of e. 


The convergents of e are 

2 3 8 11 19 87 106 193 1264 
1 ’ I’ 3’ T’ T’ 32’ ‘39’’ 71 ’ 465 ’ ' 


Using the result of the Corollary we know that the convergent lies 
1 19 1 

within-of e, — lies within-of e, and so on. The first convergent 

4 x 7 7 7 x 32 

193 

guaranteed to be accurate to 0.0001 of e is —- which lies within 

---= 0.0000303 ... of e. The previous convergent does not have 

71 x 465 f b 39 

the required accuracy because, using the left-hand inequality in the 
106 

Corollary, the error in is at least 


1 


2 x 39 x 71 


0.00018. 


The values to six decimal places 
are as follows, 
e = 2.718282 
^ = 2.718310 

^ = 2.717949 
39 


Problem 4.2 _ 

The ICF of \/2 is [1, (2)]. Determine the first convergent which is accurate 
to within 0.0001 of the true value of \/2- 


193 

We now appreciate that the convergent —— is a very good approximation 


to e, but as yet we have made little progress towards establishing the sense 
in which it is ‘best’. We shall put this right immediately by showing that of 


all rational numbers with denominator not exceeding 71, the convergent 


193 

TT 


is the best approximation. The general result from which we can draw this 
as a conclusion follows quickly from the slightly stronger result of 
Theorem 4.2 below. 


For completeness we have given proofs of all the remaining results in this 
section but you should regard these proofs as optional reading. 


Theorem 4.2 

Let — be the nth convergent of the irrational number x, and let °r be 

q n b 

any rational number with 0 < b < q n +i- Then, for n > 1, 

\q n x-p n \ < \bx-a\, 

the equality holding only if a = p n and b = q n . 


Remember — is in lowest terms. 
b 
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Proof of Theorem 4.2 

Let - be a rational number (in its lowest terms) with b < q n +i, and suppose 

it is not the case that a — p n and b = q n , when the required equality 
certainly holds. 

Consider the system of equations (with r, s real for the moment): 
a = rp n + sp n+x ; 
b = rq n + sq n+x . 

Multiplying the first equation by q n and the second by p n and subtracting 
leads to 


aq n ~ bp n — s{p n+x q n - p n q n +i) 

which, after applying Theorem 1.3 property (a) to the bracketed term on the 
right-hand side, gives 

s = (—l) n+1 (ag n - bp n ). 

Similarly 

r = (—l) n+1 (6p„ + i - aq n+x ). 


Note that r and s are therefore integers, and neither of them is 0 because 
the assumption that r = 0 leads to aq n + x = bp n + x , whereupon Euclid’s 
Lemma forces q n+x to divide b, since gcd(p n+1 ,q n+x ) = 1, and this 

CL X) 

contradicts 0 < b < q„+ x . Also if s — 0 we have — = —, which is the case of 

b q n 

equality which has been excluded. 


Now since 0 < b = rq n + sq n+x < q n +i, the non-zero integers r and s must 

have opposite sign. Moreover x — — and x — ^ w+l have opposite signs 

Qn Qn +1 

because the convergents are alternatively greater than and then less than 
the irrational limit x. It follows that r(q n x - p n ) and s(q n+x x - p n + 1 ) have 
the same sign. (We need this fact so that we can replace the absolute value 
of their sum by the sum of their absolute values.) 

I bx - a\ = | (rq n + sq n+1 )x - ( rp n + sp n+ i)| 

= | r(q n x - p n ) + s(q n+x x - p„+i)| 

= \r\\q n X-p n \ + Nkn+lZ ~Pn+l\ 

> \A\<lnX-Pn\ > \q n X-p n \ 

as claimed. ■ 


As a consequence of this result we can now spell out the sense in which the 
convergents are best approximations to the irrational number. 


Corollary Convergents are best approximations 

Pi P‘2 

Let —, —, ... be the convergents of the irrational number x. Then, 

Qi <72 

for any rational number — with 1 < b < q n 


Pn 
X - 

< 

a 

x — — 

<7n 


b 


the equality holding only if a = p n and b = q n . 


Remember y is in lowest terms. 
b 


This means for example, that as ——— is a convergent of e there is no better 
a 465 

rational approximation — to e with 0 < b < 465. 
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Proof of the Corollary 

Assume, to the contrary, that there is a rational number ^ with 1 < b < q n 
and 


a 


Pn 

x — - 

< 

X - 

b 

Qn 


Then 

\q n X~Pn\ = Qn 
contradicting the result of the theorem. 



Pn 


a 

> b 

a 

Qn 

x - 

Qn 

Qn 

x s 

x -i 


= \bx — a\, 


4.3 Rational approximations to 7r 

Having seen that there are patterns for the ICFs of all quadratic irrationals 
and also for the irrational e, we would be lacking curiosity if we did not 
enquire about the ICF of the most famous of all irrationals, ir. 
Disappointingly, it transpires that there is no known pattern to the sequence 
of partial quotients of tt. We can see how the ICF begins by working with 
some good known approximation to 7 r. For example, take 7 r = 3.141592654. 
With the aid of a calculator, the continued fraction algorithm yields the 
following. 


7 r = 3 + ( 7 r- 3)=3 + 0.141592654; 


1 

0.141592654 ' 
1 

0.0625133059 

1 

0.99659454 ~ 
1 

0.003417099 ' 


= 7.062513305; 
= 15.99659454; 
1.003417099; 

= 292.6458622; 


aj = 3 
a 2 = 7 
a 3 - 15 

«4 = 1 

tt 5 — 292 


Calculations like these are very 
sensitive to rounding so your 
calculator may well produce 
different results from ours. 


This, of course, is only giving a crude idea of how the ICF of 7 T begins. Here 
we are not calculating the ICF of tt itself but rather the continued fraction of 
the rational approximation to nine decimal places with which we began, and 
as we are not calculating the reciprocals exactly, rounding errors are 
creeping in. For example, if we were to use a calculator which works only to 
eight decimal places the ensuing value of a 5 would turn out to be 
significantly different from 292. But the partial quotients we have obtained 
thus far are correct for 7r. 


tt= [3,7,15,1,292,...]. 


The resulting sequence of convergents begins 
3 22 333 355 103 993 
T’T’ 106’ 113’ 33102 ' 

The search for good rational approximations to tt has a long history. With 
the aid of today’s computers the value of 7r is now known correct to millions 
of decimal places. The first known rational approximations can be traced 
back to Archimedes (287-212 BC). Archimedes utilized the fact that the 
circumference of the circle lies between the total edge length of any inscribed 
polygon and that of any circumscribed polygon. By using a polygon of 96 
sides he discovered lower and upper bounds given by 


223 22 

—- < 7T < —. 

71 7 

22 — 

The value —, known as the Archimedean value of tt, is commonly used as 

an approximation to 7 r, being accurate to within 0.001264. The lower bound, 
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223 . 

—is not a convergent of 7 r but turns out to be more accurate than 


22 

T’ 


the error being approximately 0.0007476. The next convergent of n, -, 

106 

must be more accurate still because we know it is the best approximation of 

333 

all rational numbers with denominator not exceeding 106; the error in —— is 
, 106 
approximately 0.000083. 


However the next convergent, ——, is remarkably accurate. (This value for 7 r 

llu 

seems to have been first known to a Chinese mathematician Tsu Chung-Chih 
(430-501 AD) but was not discovered in Europe until late in the sixteenth 
century.) Because of the large partial quotient 292 there is a huge jump in 


the values of numerator and denominator following the convergent 
the corollary to Theorem 4.1 gives the bound on its error: 


355 

113’ 


and 


1 _ 3 

113 x 33102 < 10^' 

355 

Our theory tells us that —— is a better approximation than any other 

-L J-O 

rational with denominator not exceeding 113. In fact this claim is not very 
flattering, as the following example will show. 


7 T — 


355 

H3 


Example 4.2 

355 

Show that —— is the best of all rational approximations to 7 r with 

J. J.O 

denominator not exceeding 10000. 
r 

Suppose that - is a rational number with 114 < s < 10000. Then 
s 


r 355 
s“Tl 3 = 
But we also have 
3551 


|113r - 355s| 
113s 


> 


1 


> 


113s ~ 113 x 10 000 


7 T — 


113 


< 


1 


355 

Hence —— lies within 

llu 


113 x 33102 
1 


113 x 33102 


of 7 r but is further than 


1 


113 x 10000 


from -. The triangle inequality gives 


r 

- 7 r 

> 

r 355 


355 

7 T- 

s 


s ~~ 113 


113 


> 


1 


113 x 10000 
23102 


113 x 10000 x 33102 


113 x 33 102 

> 


1 


113 x 33102 
355 


Thus - is not as good an approximation to 7r as is 


The triangle inequality for real 
numbers says 

|a — 6| < |a — a:| + |a: — 6|. 
The form used here is 

|a — x\ > |a — f>| — |a: — 6|. 


Problem 4.3 ___ 

Let x = [0,1,2,3,4,5,6,...]. Determine the first convergent of x which is 

within 10 -5 of x. Show that of all rational numbers -, where the positive 

s 

denominator s does not exceed 100, the most accurate approximation 
30 

to a: is —. 
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We have seen that the convergents to the irrational number x satisfy 

x — — < —-— and, as the denominator sequence { q n } is increasing, 

Qn QnQn+1 
'P 1 

x -— < -5-. A question we might ask concerns whether or not this 

q n qi 

property characterizes convergents; that is, if a rational number - has the 

property x — 7 < i , is 7 necessarily a convergent of x? It is not difficult 
6 6 2 6 

3 1 3 

to show that the answer is no. For example, \/3 — - < - and yet - is not 

a convergent of \/3. Nevertheless there is an inequality of this kind which 
does characterize convergents. 


Theorem 4.3 Which rationals are convergents? 

If the rational number 7 satisfies 
0 

a 1 
X ~b < 262 


then 7 is a convergent of x. 
b 


Proof of Theorem 4.3 

Suppose to the contrary that 7 is a rational number (in lowest terms) 

b Pn f 

satisfying the given inequality but which is not one of the convergents — of 

Qn 

x. Let r be the unique integer for which q r <b < q r + 1 . Then, by 
Theorem 4.2, 


| q r x - p r \ < \bx — a\ = b x - ^ < 6 x — 


~'-"r 61 - 262 26■ 

Therefore 

Pr 1 . . , . Pr 1 

q r x - < —, which gives x -< 

q r 26 q r 2q r b 

The triangle inequality then gives 


a p r p r a 11 , x 

t-— < + x-t + (*) 

6 q r q r b 2q r b 2b z 

On the other hand 

a _ Pr _ g r « ~ Prb > J_ ^ 

6 q r q r h — q r b 

(The final inequality comes from the fact that q r a — p r b is an integer and is 

CL p r 

non-zero, for otherwise 7 would be equal to the convergent —, which we 
6 qr 

have assumed not to be the case.) Combining (*) and (**) we have 

J_ _L _L 

q r b < 2q r b + 26 2 ’ 

and this simplifies to q r > 6, giving a contradiction. ■ 


Problem 4.4 - 

17 2 + y/U 

Without calculating the ICF show that — is a convergent of---. 

16 5 
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ADDITIONAL EXERCISES 


Section 1 


Express each of the following rational numbers as an FCF. 
(b) (c) ^ (d) 3.142 


23 

(a) 54 


17 


24 


2 (a) If x > 1 and x has FCF [a \, a 2 , £ 13 ,... a n show that — has FCF 

x 

[0, a\, 02 , a, 3 ,... a n \. What is the corresponding result when 


0 < x < 1 ? 

(b) Suppose that -1 < x < 0 and that x has FCF [-1, a\, a 2 ,... a n ]. 

Show that [ai, « 2 ,... a n ] has value-. 

1 + x 


3 Evaluate the sequence of convergents for each of the following FCFs 
and, in each case, verify that Theorem 1.4 holds. 

(a) [2,2,4,2,4,2] (b) [1,2,3,4,5, 6 , 7] 


4 Use continued fractions to find the general solution of each of the 
following linear Diophantine equations. 

(a) 17x + 91y = 1 (b) 16z + 21% = 17 


Section 2 


1 


Evaluate the following periodic continued fractions. 

(a) [(5,1,3)] (b) [0,(5,1,3)] (c) [1,2, (5,1,3)] 


2 


e 2/fc _ 1 

Euler showed that the ICF of the irrational number —- is 

e 2/k + i 


[0, k, 3k, 5 k, 7k, 9k ,...]. 


Use this information to find values correct to four decimal places for the 
following. 


(a) 


e — 1 
e + 1 


(b) 


e 2 — 1 
e 2 + 1 


Section 3 


Determine the ICFs of the following irrational numbers. 


(a) Vn (b) 


1 + ^13 


(c) 


4 + \/2l 


2 v ~' 5 

If 0 < k < 2n and k divides 2 n, show that 

\Jn 2 + k = n,(^-,2n 

L \ k 

Hence find ICFs for (a) \/l2, (b) \/24. 

Determine the ICF of - n, for a general n > 2. Use this to write 
down the ICFs of s/2, \/6 and \/l2. 
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Section 4 

1 The irrational number x has ICF [0,2,4,6,8,10,...]. Find the best 

rational approximation — to x with: 

b 

(a) b < 100; (b) b < 1000. 


Theorem 4.3 tells us that any rational number — satisfying 

b 


a 

X ~b 


< 


2 ft 2 ’ 


is a convergent of x. It is not true that every convergent — of x satisfies 
the inequality. 

However, prove that of two consecutive convergents — and ^ n+1 of x, 

Qn Qn +1 

at least one of them does satisfy this inequality. Hint: Assume the 
contrary and use the fact that 


Pn+1 

Pn 


_ Pn+1 

+ 

Pn 

Qn+l 

Qn 


x - 

Qn+l 

X - 

Qn 


3 A theorem of Hurwitz (1891) shows that for each irrational number x 

there are infinitely many rationals — satisfying 

b 

a 1 

x — — < _ . 

b y/h b 2 

(a) Why must any rational number ^ satisfying the above inequality 
be a convergent of xl 

(b) Find two convergents of \f?> which satisfy this inequality, and two 
which do not. 


Challenge Problems 

1 In applying the Continued Fraction Algorithm to determine the ICF of 
y/n we obtain the sequence of equations 

y/n = xi = [ai,i 2 ] = [01,02,2:3] = [o 1} a 2 ,a 3 , ar 4 ] = •••, 

where, for all k > 1, a*, = intfyfy and the step from one stage to the 
next is accomplished by 

1 

Xk — Ofc T . 

%k-\- 1 

Prove that for each k > 1, 

rk + Vn 

Xk = -, 

Sk 

where {r^} and {sfc} are sequences of integers given by 

n — rl +1 

n =0, Si = 1, nt +1 - OfeSfc - r k , S k+ 1 =-. 

Sk 

Use this sequence to determine the ICF of y/l4. 

2 Let x — [oi, 02,03,04, ...]. Find the ICF of — x. 
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3 Let x = [cti, a 2 , a 3 , a 4 ,...]. Show that provided a 2 is even, the ICF of 2x 
satisfies 

2[ffll) ^21 * * •] [2nj, 2^2j 2[fl3, CI4, tt5,. . .] ], (*) 

and deduce that provided all the even numbered partial quotients are 
even, 

2[fll, 0 2 , @ 3 > ® 4 j • ■ •] [^Q-l, 2®21 2fl3, 2U4, 2 d 5 , 2^61 ■ ■ • 

Find a formula corresponding to (*) for the case of a 2 being odd. 


4 Let = 
P3 


0 1 

1 a*. 


and M, 


_ f Qn Pn \ 

\ 9n+l Pn +1 / 


where ?i, 

9i 92 


—, ... are the convergents of [ai, a 2 , 03 , 04 ,.. .1. Show the following. 
93 

(a) M n = A„ + iM„_i, for n > 1 

(b) AF n — j4 n 3_i^4. n2 4 n _4... A 2 Aj 

(c) By considering the determinant of M n , deduce the formula of 
Theorem 1.3 property(a): 

PnQn—1 ~ QnPn—l = (-1)”, for Tl > 2. 

(d) By considering M n and its transpose, deduce the formulae 

Pn 


Pn -1 


\ a m Q-n—1 > &n— 2j •••> &l] 


and 


Qn 

Qn— 1 


— [fl n , fl n _i, tt n _2, . • • , 0-3, ^ 2 ]• 


The transpose of 


SOLUTIONS TO THE PROBLEMS 

Solution 1.1 

(a) RM-« + -\-»+5-5 

+ 4 

(b) [0,2,2,2] = 0 + L_= 0 +-^-=0+-| = ^ 

2 H-2+| 12 12 

2+ 2 

Solution 1.2 

(a) 11 = 1 x 7 + 4 

7=1x 4 + 3 

4 = 1x 3 + 1 

3 = 3x1 

Soy = [1,1,1,3], 

(b) 43 = 2 x 19 + 5 
19 = 3x5 + 4 

5 = 1x 4+1 

4 = 4x1 

So ^ = [2,3,1,4], 
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Solutions to the Problems 


Solution 1.3 

(a) 65 = 2 x 23 + 19 
23 = 1 x 19 + 4 
19 = 4 x 4 + 3 

4 = 1x 3+1 
3 = 3x1 
65 

So —= [2,1,4,1,3]. 

(b) 23 = 0 x 53 + 23 
53 = 2 x 23 + 7 
23 = 3 x 7 + 2 

7 = 3 x 2 + 1 
2 = 2x1 
23 

So — = [0,2,3,3,2], 

(c) -94 = -8 x 13 + 10 

13 = 1 x 10 + 3 
10 = 3 x 3 + 1 
3 = 3x1 
-94 

So —= [-8,1,3,3]. 

Solution 1.4 

(a) [1,2,3] = [1,2,2,1] 

(b) [-4,3,1,6,1] = [-4,3,1,7] 

(c) [0] = [—1,1] 

Solution 1.5 

(a) The numerators of [2,3,4,1,2] are 

Pi = 2; P 2 = 2x 3 + 1 = 7; p 3 = 4 x 7 + 2 = 30; 

Pa — lx 30 + 7 = 37; p 5 = 2 x 37 + 30 = 104. 

The denominators are 

9i = l; q 2 = 3; q 3 = 4 x 3 + 1 = 13; 

94 = 1 x 13 + 3 = 16; q 5 = 2 x 16 + 13 = 45. 

The convergents are therefore 
2 7 30 37 , 104 

1’ 3’ 13’ 16 and ~45~' 

(b) The numerators of [—2,3,1,6] are 

Pi = -2; p 2 = (-2) x 3 + 1 = -5; 
p 3 = 1 x (-5) + (-2) = -7; p 4 = 6 x (-7) + (-5) = -47. 
The denominators are 

q\ = 1; 92 = 3; 93 = 1 x 3 + 1 = 4; 94 = 6 x 4 + 3 = 27. 

The convergents are therefore 

-2 -5 -7 J -47 
——, , —— and ——. 


Note that the integer part of — — 
is —8 (not —7). 
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Solutions to the Problems 


Solution 1.6 


k 

1 

2 

3 

4 

5 

Pk 

-4 

-3 

-19 

-22 

-151 

dk 

-4 

1 

5 

1 

6 

Qk 

1 

1 

6 

7 

48 


The convergents are 


(b) 


-4 -3 -19 -22 -151 

1 ’ 1 ’ "IT’ T and ~48~ 


k 

l 

2 

3 

4 

5 

6 

7 

Pk 

l 

2 

5 

7 

19 

26 

71 

a k 

l 

1 

2 

1 

2 

1 

2 

qk 

l 

1 

3 

4 

11 

15 

41 


The convergents are 


1 2 5 7 19 26 71 

1’ 1’ 3’ 4’ Tl’ 15 and 41- 


Solution 1.7 


(a) The convergents of [2,3,4,1,2] were found to be 
2 7 30 37 J 104 
1’ 3’ 13’ 16 and 

The FCF [2,3,4,1, 2,2] has the same first five convergents plus a sixth, 
which is its value. The latter is determined from the previous 
convergents as 

2 x 104 + 37 245 

2 x 45 + 16 ~~ 106" 


(b) Similarly [2,3,4,1,2,100] has the same first five convergents as 
[2,3,4,1,2] and has final convergent 

100 x 104 + 37 10437 

100 x 45+ 16 ~~ 4516 ' 


The addition of a further, even larger, partial quotient has not changed 

104 

the value of the FCF by much. In fact -— 2.3111111, 

^ - 2.3113208 and - 2.3111160. 

106 4516 
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Solutions to the Problems 


Solution 1.8 

From the Euclidean algorithm 
48 = 2 x 19 + 10 
19 = 1 x 10 + 9 
10 = 1 x 9 + 1 
9 = 9x1 
48 

and so — has FCF [2,1,1,9]. Its convergents are 
2 3 5 .48 

I' I' 2 am> S' 

From the final two convergents we have 
48x2-5x19 = 1 

and therefore, multiplying through by 4 and rearranging, 

48 x 8 - 19 x 20 = 4. 

T his gives x = 8,y = — 20 as one solution of the linear Diophantine equation 
48a: + 19y = 4, and so the general solution is 

x = 8 + 19f, y = -20 - 48f, t € Z. 


Solution 1.9 

The convergents of [2,1,1,2,1,1,2] are 

2 3 5 13 18 31 ,80 

—, —, — and —. 

1’ 1’ 2 5 7 12 31 

Theorem 1.4 orders the convergents as 


5 18 80 31 13 „ 

2< 2 < Y < 5I < T2 < y <3 ' 


which is readily confirmed by evaluating these as decimals. Our results are 
recorded to a maximum of three decimal places. 


2 < 2.5 < 2.571 < 2.581 < 2.583 < 2.6 < 3. 


Solution 1.10 

For any FCF all partial quotients (with the possible exception of ai) are 
positive. So q\ = 1, 92 = a 2 and q% = 0302 + 1 satisfy 1 = 91 < 92 < 93 - 

Heading for an induction proof, suppose that % > qk -1 > 1 for some k > 3. 
Then 


To avoid approximate division note 
that, for example, 5 X 31 < 12 X 13 
31 13 

from which — < — follows. Each 
12 o 

of the inequalities can be tested in 
this way by cross-multiplication. 


qk +1 = o,k+iqk + Qk -1 > Qk + Qk -1 > Qk + 1 > Qk, 


and the induction step is established and the proof is complete. 


Finally, since <? 3 > 2, the above inequality qk+i > Qk + Qk- 1 shows that, 
from qq onwards, the qkS increase in steps of at least 2 and as qq > 3, qk > k, 
for k > 5 follows. (This could be proved slightly more formally by using 
Mathematical Induction.) 


Solution 2.1 

The convergents of the ICF [1,1,1,!,...] begin 


1 2 3 5 8 13 21 
I’ 1’ 2’ 3’ 5’ IT’ 13' 


The error in C 5 = - is bounded above by --- = 0.025. 

5 5x8 

13 1 

The error in Cq = — is bounded above by --— = 0.009615 (approx). 

8 o X lo 

The actual error in C 5 is approximately 1.61803 — 1.6 = 0.01803. 

The actual error in C$ is approximately 1.61803 — 1.625 = —0.00697. 


In fact, the ‘smallest’ sequence of 
denominators occurs when every 
partial quotient (after a 1 ) takes its 
minimum value of 1 , and then the 
sequence of denominators is 1 , 1 , 2 , 
3, 5, 8 , ..., which is the Fibonacci 
sequence. 
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Solutions to the Problems 


Solution 2.2 

If we put x = [3,1,3,3,1,3,3,1,3,...] then, as in Example 2.1, 
x = [3,1,3, x\. The convergents of [3,1,3, x\ are 

3 4 15 , 15a:+ 4 

V V T and toTT' 

and, as the final convergent is x itself, 

15a; + 4 9 

x = — -—, that is 4x — 14x — 4 = 0. 

Ax + 1 

As x > a\ = 3, it is the positive root of this quadratic equation and 
7 +^65 

SO X = ---. 


Solution 2.3 

2 

Putting x = [2,2,2,2,...] we have x = [2, x]. The convergents of [2, a;] are - 
2x + 1 . 

and -. As the second of these is equal to x itself we reach the 

x 

quadratic x 2 — 2x — 1 = 0 which has roots x = 1 ± \/2. As x > ai = 2 it is 
the positive root of this quadratic equation and so 

[2, 2, 2, 2,...] = x = 1 + y/2. 

Now notice that we can add an integer to the value of an ICF by adding it 
to the leading partial quotient ai, that is 

[ai,a 2 , a3, ■ ■ •] + n = [ai + n,a 2 ,a 3 ,...]. 

Here, adding -1 to both sides, 

[1,2,2,2,...] = V2, 

which gives the ICF of \[2. Its convergents begin: 

1 3 7 17 41 99 239 
1’ 2’ 5’ 12’ 29’ 70’ 169' 

Now \[2 lies between each pair of successive convergents. So, for example, it 

lies between — and —, whose difference in size is ---. Hence — lies 

29 70 29 x 70 29 


within 


1 


29 x 70 V 


which is less than- 


1 


1000 


of y/2. 


Solution 2.4 

(a) Let x = [(1,1,1,4)] = [1,1,1,4, x}. The convergents are 
1 2 3 14 14a; + 3 

1 ’1’2’y and ~9^+2' -x ' 

The resulting quadratic equation is 9x 2 - 12x -3 = 0 which, after 
cancelling the common factor 3, becomes 3.x 2 - 4x - 1 = 0. This has the 
2 + y/7 

positive root x = -. So 

3 

[2, <1,1,1,4)] = [2, x] = 2 + — = 2 H- 

x 2 + V7 


7 + 2y/7 = y/7x (V7 + 2) r- 
2 + y/7 2 + \/7 


(b) Let x = [(3)] = [3,x]. Then x = 3 H—, giving the quadratic 

X 

3 +x/13 


x 2, — 3x — 1 = 0 which has positive root x = 


Now 


[2,2, <3)] = [2,2, x\ 
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Solutions to the Problems 


2 5 5x + 2 

has convergents -, - and n i 1 , and therefore 


[2,2,(3)] = 


12 2x + 1 
5x + 2 


„ 3 +Vl3 

5 x - 2 - +2 _ 19 + 5VT3 


2x + 1 3 + s/13 8 + 2y/l3 ’ 

2 x---+ 1 


and this, when simplified by multiplying top and bottom by 8 — 2s/l3, 
gives 

11 + -/I3 


[2,2, (3)] 


6 


Solution 3.1 

(a) x\ = s/3 lies between 1 and 2 and so int(-\/3) = 1. Therefore 

1 


s/3 = 1 + (Vz - 1), giving m = 1 and X 2 = 

Repeating for X 2 we have 
1 


Vz-i 


Vz + 1 Vz-I . . , . 2 

—-— = H--—, giving a 2 = 1 and x 3 = 


s/3-1 


s/3-1 2 

From £ 3 : 

2 1 

—-= s/3 + 1 = 2 + (s/3 - 1), giving a 3 = 2 and x 4 = — 7 = —r- 

s/3-1 v 3 — 1 

As x 4 = x 2 we are now in a repeating cycle with a 4 = a 2 , 0.5 = a 3 , 
as = 04 , etc. So s/3 = [1, (1,2)]. 


1 


x 2 = 

x 3 = 

x 4 = 

x 5 = 


s/ 7-2 s /7 + 2 

3 s/7+1 


s/7-1 s/7 + 1 

_2_ s/7 + 1 

s/7-1 X s/7+1 

3 s/7 + 2 r- 

x — 7 =-= v 7 + 2 = x\ 


2 ); 



ai = 4, 

X2 = 

s/7 + 2 
3 

= 1 + 

y/7-1 

3 ’ 

a 2 = 1, 

x 3 = 

s/7+1 

2 

= 1 + 

s/7-1 

2 ' 

a 3 = 1, 

X4 = 

s/7 + 1 
3 

= 1 + 

“7* 

to 

0,4 — 1, 

x 5 = 


1 


s/7-2 

3 

y/7-1 

2 

y/7-1 

3 

s/7-2 


s/7-2 s/7 + 2 

Hence a 3 = 01 , - a 2 , 0,7 = a 3 , etc. 

Therefore 2 + s/7 = [(4,1,1,1)]. 

Solution 3.2 


For n > 1 the integer part of sjn 2 + 2 is n and so we write 

s/n 2 +2 = n + (s/n 2 + 2 — n), 

1 


giving a\ = n and x 2 = 


s/n 2 + 2 — n 
multiplying top and bottom by s/n 2 + 2 + n: 

s/n 2 + 2 + n 


. This expression for x 2 is simplified by 


1 

s/n 2 + 2 — n 
giving a 2 = n and x 3 = 
2 


n + 


s/n 2 + 2 


s/n 2 + 2 - n 
giving a 3 = 2n and x 4 = x 2 - 


s/n 2 + 2 - n 
s/n 2 + 2 + n = 2nd- (Vn 2 + 2 - n), 


In each of the last two calculations 
we multiplied top and bottom by 
s/3 + 1 and then pulled out the 
integer part. 
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Solutions to the Problems 


This pair of equations now cycle and so y/n 2 + 2 = [n, ( n , 2n)]. 

Hence we have: 

V3= [1,(1,2)] ; 
v/6 = [2, (2,4)]; 

^11 = [3, (3,6)]; 

V^8 = [4, (4,8)]. 

Solution 3.3 

(a) From the evidence presented in Table 3.1 you might have made some of 
the following observations: 

• the ICF of y/n is periodic, with cycle starting at 0 , 2 ; 

• the final partial quotient in the cycle is 2aj; 

• leaving aside the final term 2ai, what remains of the cycle is 
symmetric (palindromic), reading the same forwards and backwards; 

• the palindromic part sometimes has even length, sometimes odd, 
and in some cases (e.g y/2 and \/5) is non-existent. 

(b) The effect is as follows. 

Adding ai, which is equal to the integer part of y/n, gives a first partial 
quotient of 2ai, and this can be absorbed into the cycle leaving a purely 
periodic ICF. For example, 

V7 + iDt(>/7) = [(4,l,l,l)]. 

(c) Collecting all this together we might make the conjecture that the ICF 
of y/n takes the form 

y/n= [ai,(a 2 ,a 3 ,... 03,02,2(2!)]. 


Solution 4.1 

The convergents of y/5 are 
2 9 38 161 
1’ 4’ 17’ ~72’ 

The convergent ^ satisfies 


1 


< 


38 
X ~ 17 


< 


1 


17 x 72' 


2 x 17 x 72 
That is, working to five decimal places 
38 | 


0.00041 < 


x — 


17 


< 0.00082. 


Solution 4.2 

The convergents of y/2 are 

1 3 7 17 41 99 239 
1’ 2’ 5’ 12’ 29’ 70’ 169’ 
99 

The convergent — satisfies 
70 


99 

70 


< 


1 


< 


1 


70 x 169 10000 


= 0 . 0001 , 


and so has the required accuracy. It is the first convergent with this 

41 

property because the previous convergent — satisfies 


41 

29 


> 


1 


> 


1 


2 x 29 x 70 10 000 


= 0 . 0001 . 
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Solutions to the Problems 


Solution 4.3 

The convergents of x are 

0 1 2 7 30 157 972 


1’ 1’ 3’ 10’ 43’ 225’ 1393’ 


157 

The convergent —— satisfies 

ZZD 


x — 


157 

225 


< 


1 


225 x 1393 


< 0.000004, 


and so has the required accuracy. It is the first convergent with this 

30 

property because the previous one, —, satisfies 


30 

43 


> 


1 


2 x 43 x 225 


> 0.00005. 


r 30 

Let - be a rational different from — with 1 < s < 100. Then, since 
s 43 


43r — 30s is a non-zero integer, we have 


But 


x — 


30 

43 

30 

43 


1 


< 


|43r — 30s | ^ 1 ^ 

43s “ 43s - 43 x 100' 


1 


43 x 225’ 


and so the triangle inequality gives 


> 


> 


1 


125 


43 x 100 
1 

43 x 225 


> 


43 x 225 
30 


43 x 100 x 225 


x — 


43 


30 r 

and so — is a better approximation than -. 
43 s 


Solution 4.4 

17 2 + %/n 

By Theorem 4.3, a sufficient condition for — to be a convergent of--- 


is that 


2 + Vn 17 
~5 16 


< 


1 


2 x 16 2 


This can be confirmed by use of a calculator or 


2 + Vn 17 

~5 16 


| - 53 + 16VTT| 

80 

|(—53 + 16VTI)(-53 - 16y/IT)l 
80 x |(—53 - 16 a/TT)| 

7 


80 x (53 + 16 -v/IT) 


< 


< 


1 


8000 2 x 16 2 


53 2 = 2809 
11 x 16 2 = 2816 


53 + 16vTT > 53 + 16 x 3 > 100 
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SOLUTIONS TO ADDITIONAL 
EXERCISES 


Section 1 


1 (a) 23 = 0 x 54 + 23 

54 = 2 x 23 + 8 
23 = 2 x 8 + 7 
8= 1 x 7+1 
7 = 7x1 
23 

So — = [0,2,2,1,7], 

(b) 117 = 6 x 17+15 

17 = 1 x 15 + 2 

15 = 7 x 2 + 1 

2 = 2x1 

„ 117 r 
S° ~jy = [6, 1, 7,2]. 

(c) -37= (-2) x 24 + 11 

24 = 2 x 11 + 2 
11 = 5 x 2 + 1 
2 = 2x1 
—37 

So — = [-2,2,5,2]. 

(d) We treat 3.142 as 51^?. 

V ’ 1000 

3142 = 3 x 1000 + 142 

1000 = 7 x 142 + 6 

142 = 23 x 6 + 4 

6 = 1x 4 + 2 

4 = 2x2 

So 3.142 = [3,7,23,1,2], 


Notice that this works without 
reducing the rational numbers to 
lowest terms. We leave you to 
ponder why. Generally we consider 
our rational numbers to be in 
lowest terms. 


(a) If x = [ai,a 2 , 03 ,... a n ] > 1 then 

1 


1 „ 1 „ 

— — oh— — 0 + -— 

X X [GO, 0 2 1 O 3 , . . . 0„ 


[ 0 , Oj\ , 0 * 2 , CL 3 , • - - drj 


Note that aj is not the leading partial quotient in this FCF and so 
it is essential that (ti > 1, a fact which is guaranteed by x > 1. 

For 0 < x < 1 we have — > 1, and so we can apply the result just 

obtained to —. If — = [&i, b 2 , 63,... b n ] then the FCF of its 
x x 

reciprocal (noting that = x) is given by x = [0,61,6 2 ,63 ,... b n }. 

In other words, if 0 < x < 1 then its FCF will have leading partial 

quotient 0, say x = [0, oj, a 2 ,...a n ] and then - = foi,a 2 ,... a„]. 

x 

(b) If x = [-l,oi,a 2 , • - .an] then 1 +x = [0,a 1; a 2 ,... a n }. By the 
result of part (a) the reciprocal of 1 + x has FCF [01, a 2 , ... a„]. 


50 





(a) [2,2,4,2,4,2] has convergents 

2 5 22 49 218 485 

I’ 2’ ¥’ 20’ ~89~ and 198 

According to Theorem 1.4 we have the inequalities 

2 22 218 485 49 5 

T < T < "89" < 198 < 20 < 2’ 
each of which is readily verified. 

(b) [1,2,3,4,5,6,7] has convergents 

1 3 10 43 225 1393 9976 

T’ 2’ T’ 30’ 157’ 972” ““ 6961' 

The inequalities given by Theorem 1.4 are 

10 _ 225 9976 ^ 1393 ^ 43 _ 3 

1 < Y 157 < 6961 < "972" < 30 < 2’ 

each of which is readily verified. 


(a) The FCF of — is [5,2,1,5] and it has convergents 
5 11 16 , 91 

I- T' 3 and Tr 

From the final pair of convergents 3x91 — 16x17=1, which 
shows that x = —16, y = 3 is one solution of 17x + 91 y = 1. The 
general solution is therefore given by 

a: = -16+ 91 1, y = 3-17t, teZ. 


(b) The FCF of yy is [13,1,2,5] and has convergents 
13 14 41 , 219 

T' T’ T “<* w" 

From the final two convergents we have 3 x 219 — 16 x 41 = 1. 

Multiplying through by 17 gives x = (—41) x 17 = —697 and 
y = 3 x 17 = 51 as one solution of 16.x + 219 y = 17, and the general 
solution is given by 


x = -697 + 219f, y = 51 — 16f, t € ~L 
Alternatively since t = 3 gives the solution x 
x = -40 + 219f, y = 3 - 16t, t e Z. 


-40, y = 3, 


Section 2 


1 (a) Let x = [(5,1,3)]. Then x = [5,1,3, x] and this has convergents 

5 6 23 , 23x + 6 

— and-. 

1’ 1’ 4 4x + 1 

As the final convergent is equal to x we have the quadratic equation 


4x 2 - 22x - 6 = 0. 


The positive root of this equation is x 


11 + >/l45 
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(b) [0, (5,1,3)] has value —, where x is as in part (a). That is 


[0, (5,1,3)] 


and th is ca n be simplified by multiplying top and bottom by 

11 - VT45: 


[0, (5,1,3)] 


11 + -/L45’ 
iplified by m 

4(11 -x/145) 7145-11 


-24 


6 


1 Q Q T il 

(c) [1,2, (5,1,3)] = [1,2, x] and has convergents -, - and --. So, 

12 2x + 1 

substituting for x from (a), 


[1,2, (5,1,3)] = 


which can be simplified to 


37 + 3y/l45 
26 + 27145’ 

23 + 7145 


24 


(a) Taking k = 2 we have-- = [0,2,6,10,14,18,...]. The 

e t 1 

convergents begin 

0 1 6 61 860 

T’ 2’ 13’ 132’ 1861’ ' 

61 

The convergent —- is accurate to within ——-which is less 

132 132 x 1861 

than 5 x 10 -6 and guarantees correctness to four decimal places. 


(b) Taking k = 1 we have 


are 


e 2 — 1 
e 2 + 1 


[0,1,3,5,7,9,...]. The convergents 


0 1 3 16 115 1051 
1’ 1’ 4’ 21’ 151’ 1380’ ' 

115 1 

The convergent —- is accurate to within- which is less 

151 151 x 1380 

than 5 x 10 -6 and guarantees correctness to four decimal places. 


Section 3 


1 (a) *1 = 711 = 3 + (Til - 3); 

1 


711 + 3 


Tii-3 
2 

7II-3 


3 + 



oi = 3, 

x 2 = 

3 

> 

02 = 3, 

X 3 = 

3); 

to 

II 

CO 

<3 

X 4 = 


1 


Til-3 
2 

Tii-3 
Tii-3 


As X 4 = X 2 we are now in a cycle with the second and third 
equations recurring indefinitely, so Til = [3, (3,6)]. 


1 + Ti3 „ , Ti3-3 
(b) X! =---= 2 + 


(L\ — 2, 


2 2 

7i3 + 3 Ti3-3 

3 H ---; 0,2 — 3 , 


X2 


x 3 = 


Ti3-3 
2 


Ti3-3 2 2 • — I"’ v /i3_ 3 

As x 3 = x 2 we are immediately into a cycle and, after the initial 

partial quotient of 2, the partial quotient 3 recurs indefinitely. 

Hence = [2, (3)]. 
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ai = 1 


. 4+V21 , , V21-1 

(°) *. = — 5 — = 1 + — 5 —: 

5 _ y/21 + l _ y/21 — 3 

4 + 4 5 

4 ^ + 3 V^l-3 

- = —- =9 4- —-• 

V21-3 3 3 ’ 

3 v/21 + 3 , y/21-1 

-7=- = —; = 1 + - ;-; 


= 1 + 


y/21 — 1 
4 

4 

v/21+3 

y/21 — 3 

3 

3 

v/2l+3 

v/21-3 

4 

4 

v /21 + 1 

\/ 2 T — 1 

5 

5 

v^I-4 

= \/2l + 4 

1 

v^l + 4 


a 2 = 1 , x 3 = 

a 3 = 2 , X 4 = 

(z 4 — 1, X5 — 

d 3 — 1 , Xq — 


5 

y/21- 1 

4 

V^l-3 

3 

v^I-3 

4 

v/2l- 1 

5 

y/21-4 


- 7f 5_^=V2T + 4 = 8 + (V2l-4 ); «=8, x 7 = ^ 

1 y/21 + 4 V2T-1 , 5 

V2T-4 5 + 5 ’ “ 7 *’ X8 V21-1 

Now ^8 = x 2 and we are in a cycle, with the sequence 1,2,1,1,8,1 
of partial quotients recurring. So 

= [1, (1,2,1,1,8,1)]. 

But notice that this gives 

= [(1)1)2,1,1,8)]. 

In fact, although it was not immediately obvious, x-j = x\ and we 
were in a cycle at the point of calculating Xj. 

2 Since 0 < k < 2n we have n 2 < n 2 + k < n 2 + 2n + 1 = (n + l) 2 . 
Therefore the integer part of y/n 2 + k is n, and so 

/ ii " 7* / / n 7 s 1 


'n 2 + k = n + (v n 2 + k — n); cq = n, x 2 = 

1 \/n 2 + k + n 

+i 2 + k — n k 


i 2 + k — n ’ 


Vn 2 + k — n k 

The integer part of y/n 2 + k + n is 2 n, and as k divides 2 n we have 

1 2 n y/n 2 + k — n 2 n k 

y/n 2 + k — n k k k y/n 2 + k — n' 

k / „ . _ , r ~ k 


/n 2 + k — n 


= y/n 2 + k + n = 2 n + (y/n 2 + k — n); 03 = 2 n, 


y/n 2 + k — 1 
2 n 

As x 4 = x 2 we are in a cycle with the pair of partial quotients —, 2 n 

k 

repeating indefinitely. So 


y/n 2 + k = n, ^ —, 2 nyJ. 

(a) Take n = 3 and k = 3. Then y/l2 = [3, (2,6}]. 

(b) Take n = 4 and k = 8. Then y/24 — [4, (1,8)]. 




3 For n > 2 we have 


(n - l ) 2 = n 2 - 2 n + 1 < n 2 - n < n 2 
and so \Jn ' 2 — n lies strictly between n — 1 and n. Therefore 

X\ = \Jv? —n = n — 1 + (y/ n 2 — n — (n — 1)), giving a\ = n - 1. 
Continuing the Continued Fraction Algorit hm - 

__1__ Vn 2 — n + (n — 1 ) \/n 2 — n + (n — 1) 

\/n 2 — n — (n — 1) (n 2 — n) — (n — l) 2 n — 1 

A s the in teger part of \Jn 2 — n is n — 1 the integer part of 

Vn 2 - n + n - 1 

--- is 2 and so 

n — 1 


x 2 = 2 + 


Vn 2 — n — (n — 1) 
71—1 ’ 


giving a 2 = 2 . 


At the next stage we look at 
ri—l 


X3 


yjn 2 — n — (n — 1) 

As this has integer part 2(n — 1), 


s/n 2 — n + (n — 1). 


x 3 = 2 (n — 1) + (y/n 2 - n — (n - 1)), giving a 3 = 2 (t7 — 1). 

And there is good news. The fractional part, \fn 2 - n - (n - 1), is the 
same as the fractional part of Xi, so we are in a cycle and can conclude 
that 


y/n 2 — n = [n - 1, (2, 2 n - 2)]. 


For example: 


n = 2 gives \[2 = [1, (2,2)] = [1, (2)]; 
ti — 3 gives V 6 = [2, (2,4)]; 

77 = 4 gives Vl2 = [3, (2,6)]. 


Section 4 

1 (a) The convergents of x = [0, 2,4,6,8,10,...] are 

0 1 4 25 204 2065 
1’ 2’ 9’ 56’ 457’ 4626’ 


25 

The convergent — satisfies 
56 


25 

56 


< 


1 


and is known to 


56 x 457 

be the best rational approximation amongst rationals with 

denominator not exceeding 56. Now suppose that y is any rational 

b 

25 

with b < 100 and different from —. Then 

56 


25 

56 


|56a — 256| 1 1 

566 “ 566 “ 56 x 100' 


It follows (from the triangle inequality) that 


a 

- - X 

> 

a 

25 


25 

x - 

b 


b 

” 56 


56 


1 1 357 ^ 1 

> 56 x 100 56 x 457 “ 56 x 100 x 457 > 56 x 457’ 

25 cl 

and therefore — is a better approximation than -. 
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(b) In the same vein, the convergent ~~ satisfies 


| 4571 457 x 4626 

So if - is a rational number different from and with 6 < 1000 
b 457 

then 


a 204 
6 _ 457 


457a - 2046 
4576 


> _L >-i_— 

“ 4576 - 457 x 1000’ 


a 

b~ x > 


457 x 1000 457 x 4626 457 x 4626 x 1000 457 x 4626' 


Therefore, - is the best approximation amongst rationals with 

457 

denominator not exceeding 1000. 

Since the convergents of x are alternately greater than and smaller 
than x, 

Pn +1 I I Pn I IPra+l Pn I \Pn+lQn PnQn+1 \ 


Pn 

+ X - 

Qn 


| <7n+l Qn 


QnQn+1 


QnQn +1 


So if neither convergent satisfies the given property 


QnQn+l 




T his can be rearranged as ( q n +\ — q n ) 2 < 0. As a square cannot be 
negative and q n +\ > q n except for the possibility of q\ = r /2 = 1, this 
gives the required contradiction for all except this exceptional case. 
However, when qi = q<z = 1 the first two convergents of x will be int(x) 
and int(x) + 1 and the result makes the modest, and true, claim that 

one of these lies within - of the true value. 


(a) Since y/b > 2, x — r < — < ttpt , and Theorem 4.3 guarantees 

w 6 y/5b 2 26 2 

that — is a convergent of x. 
b 

(b) \/3 = [1, (1,2)] and has convergents 

1 2 5 7 19 26 71 

I’ 1’ 3’ 4’ TT’ 15’ 41’ "' ' 

It can be checked by calculator that I >/3 — l| > —p, |\/3 — 21 < 

V5 v 5 

r- 5 1 , /s 7 1 , 

v/3-> —and \/3 —- < - 7 =, so the convergents - and - 

3 9^5 4 l6\/5 1 3 

2 7 

fail to satisfy the inequality while the convergents - and - do satisfy it. 
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